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S u mm  a r y 

This  is  the  final  report  describing  the  work  performed  on 
Contract  No.  F44620-74-C  - 0066  during  the  period  April  1,  1974  to 

April  1,  1975.  The  purpose  of  this  effort  is  to: 


1. 

Establish  a classification  technique  for  mixed 
seismic  events, 

2. 

Establish  a computer  model(s)  for  these  classi- 
fication techniques. 

3. 

Evaluate  the  azimuthal  and  polarization  filters 
and  array  processing  methods  as  a combined 
data  processing  technique  and  its  usefulness 
as  a discriminate  criterion. 

Seismic  verification  of  underground  nuclear  explosions  requires 
the  classification  of  seismic  events  from  observations  made  at 
teles eismic  distances.  For  single  events,  considerable  success  has 
been  achieved  by  examining  the  characteristics  of  the  body  and  sur- 
face waves  and  applying  such  single  event  classifiers  as  the  M^m^ 
ratio.  For  mixed  events,  the  key  problem  is  to  free  the  surface 
waves  of  the  event  from  the  overlapping  events  and  interferences  to 
permit  single  event  classifiers  to  be  used. 

The  decomposition  of  mixed  events  into  the  multipath  and 
interference  components  has  been  done,  in  the  past,  by  groups  unde] 
Lynn  Sykes  and  S.  S.  Alexander,  using  such  techniques  as  the  Polari- 
zation and  Azimuthal  Filters.  For  polarization-type  filters,  the 
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basis  of  the  decomposition  is  a phase  difference  sort,  such  as 
the  phase  difference  between  the  radial  horizontal  and  vertical  com- 
ponents of  a Rayleigh  wave.  For  azimuthal -type  filters,  the  basis 
of  the  decomposition  is  a frequency  domain  sort  based  on  the  ampli- 
tude ratio  of  the  horizontal  seismometer  signals.  Azimuthal  filtering 
can  also  be  based  on  array  techniques  such  as  the  high  resolution 
f-k  methods,  but  these  are  not  considered  here. 

A performance  evaluation  of  the  Polarization  and  Azimuthal 
Filters  show's  that  they  can  have  only  limited  success  in  a mixed 
event  situation  because: 


1. 


They  are  essentially  sorting  and  not  filtering 
techniques. 


2. 


The  multipath  components  of  the  mixed  events 
must  be  non-overlapping  in  frequency-time  space. 


To  avoid  these  limitations,  the  Ellipticity  Filter  was  proposed. 
The  Ellipticity  Filter,  because  the  ellipticity  of  the  Rayleigh  wave 
components  is  used  as  a constraint,  is  a structural  approach  which 
permits  quantified  relations  to  be  written  between  the  amplitudes, 
phases,  and  azimuths  of  the  multipath  signals.  Perturbations  due 
to  unstructured  components  such  as  noise  can  be  minimized  by 
means  of  auxiliary  relationships.  Moreover,  matched  filter  process- 
ing can  be  included  by  interrelating  the  components  at  different 
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frequencies.  The  overall  result  is  a unified  processing  technique 
which  should  prove  much  superior  to  the  fragmented  approaches 
previously  used. 

The  superiority  of  the  Ellipticity  Filter,  however,  requires 
accurate  calibration  of  the  seismological  stations  in  both  amplitude 
and  phase  shift  as  a function  of  frequency.  Until  recently,  the 
presence  of  substantial  drift  in  the  frequency  response  of  the  instru- 
mentation prevented  attainment  of  accuracies  much  better  than  ten 
percent.  A considerable  effort  was,  therefore,  necessarily  directed 
toward  the  attainment  of  accurate  calibration  in  the  presence  of  such 
drift.  This  effort  was  not  successful  either  because  the  sites  were 
not  calibrated  frequently  enough  for  the  complex  transfer  function  to 
be  determined  or  because  the  drift  was  unpredictable. 

Since  this  situation  has  now  been  alleviated  by  the  development 
and  cite  emplacement  of  better  electronics,  it  is  recommended  that 
the  Ellipticity  Filter  technqiue  be  tested  using  current  instrumentation. 
The  results  should  be  not  only  a good  method  for  coping  with  the 
mixed  event  situation,  but  lead  to  a better  seismological  understanding 
of  earth  physics  through  the  more  accurate  measurement  of  ellipticity. 
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SECTION  I 


PERFORMANCE  EVALUATION  OF 


AZIMUTHAL  AND  POLARIZATION  FILTERS 
FOR  SEISMIC  APPLICATIONS 


The  decomp  jsition  of  mixed  events  into  the  multipath 
components  has  been  done,  in  the  past,  using  such  tech- 
niques as  the  Polarization  and  Azimuthal  Filters.  The 
performance  of  these  filters  is  considered  here  in 
detail.  The  evaluation  shows  that  Polarization  and 
x zimurhal  Filters  are  limited  in  a mixed  event  situation 
because  they  do  not  fully  exploit  the  structure  of  the 
horizontal  and  vertical  components  of  the  multipath. 

The  result  is  that  these  techniques  can  be  successful 
only  if  the  events  are  non- overlapping  in  frequency- 
time space.  For  this  reason,  the  Ellipticity  Filter 
is  superior. 
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1.  0 


Introduction 


A mixed  event  situation  occurs  when  the  surface  seismic 


waves  of  a target  event  arrive  at  a seismological  station  at  the 
same  time  as  the  seismic  waves  of  other  events,  such  as  the  coda 


of  an  earthquake.  In  some  mixed  situations,  the  arrivals  are  not 
simultaneous  at  all  stations  and  the  events  can  be  resolved.  In  other 


situations,  however,  the  mixed  event  itself  must  be  resolved  using 
the  composite  event  data.  The  purpose  of  this  decomposition  of  the 
composite  event  is  to  determine  such  features  as  the  magnitude  of 
the  target  event  so  that  classifiers  such  as  the  Mg  : m^  ratio  can  be 
used  effectively. 

The  decomposition  of  a mixed  event  can  be  easily  done  when 
the  component  events  fortuitously  do  not  occupy  the  same  frequency 
bin  during  the  same  time  interval.  In  such  instances,  the  basic 
events  can  be  separated  from  each  other  when: 


1. 


there  are  physical  differences  between  the  events 
that  can  be  used  as  the  basis  for  the  sort. 


2. 


the  background  noise  level  is  low  enough  so  that  the 
sorting  process  is  not  excessively  interfered  with. 


Physical  differences  that  can  be  used  as  a basis  for  a sorting 
procedure  are  polarization  and  azimuth.  These  are  encompassed  in 
such  techniques  as  Choy  and  McCamy's  Polarization  and  Azimuthal 


, 
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Filters^*  Similar  techniques  are  encompassed  in  the  REMODE 

3,  4 

system  of  Sax  and  Mims  as  well  as  in  non- seismic  systems  such 


as  those  associated  with  radio  and  underwater  sound  systems. 

It  is  to  be  emphasized  that,  in  order  for  these  techniques  to 
be  successful.,  the  events  which  are  mixed  must  not  occupy  the  same 
frequency  bln  at  the  same  time.  In  addition,  the  background  noise 
level  must  be  low  to  permit  the  necessary  non-linear  processing. 

The  magnitude  of  these  problems  can  be  roughly  quantified 
using  general  principles.  At  the  decision  point  of  the  sort,  the  non- 
linear sorting  op<  ‘•ation,  because  it  is  a non-linear  process,  requires 
a signal-to-noise  ratio  threshold  of  10  to  20  dB  for  the  sorting  to  be 
successful.  That  is,  the  ratio  of  the  desired  event  power  level  to 
the  background  noise  power  level  on  a per  frequency  bin  basis  should 
exceed  about  10  dB  for  any  reasonable  performance,  and  20  dB  for 
excellent  performance.  Below  10  dB,  one  usually  has  the  problem 
of  deciding  if  the  performance  is  good  enough,  but  above  20  dB  most 
non-linear  operations  are  relatively  clean.  Now,  if  there  are  about 
a half  dozen  independent  frequency  bins  over  the  frequency  range  of 
interest,  this  implies  that  the  desired  signal  to  background  noise 
level  should  exceed  about  4 dB  for  the  sort  of  reasonable  perfor- 
mance and  12  dB  for  excellent  performance.  These  results  will  be 
more  fully  quantified  later. 
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The  requirement  of  a good  signal  to  background  noise  level 

can  be  alleviated,  to  some  extent,  by  using  the  long  period,  high-gain 

seismometer  stations  developed  by  Larmont -Doherty  Geological  Obser- 
5 

vatory  . These  long-period,  high-gain  seismometer  stations  have 

their  peak  magnification  centered  about  the  earth  noise  minimum  near 

40  seconds  and  permit  observation  of  very  long-period  surface  waves 
f 6,7,8 

from  about  25  to  60  seconds  . As  few  as  three  of  these  stations 

are  capable  of  detection,  by  direct  instpection  alone,  roughly  60  per- 
cent of  all  events  .m  a world-wide  basis  whose  magnitudes  are 

9 

greater  than  m^  = 4.  0 . With  such  a reduced  noise  background,  the 
remaining  difficulty  is  primarily  the  mixed  event  problem  or  the 
detection  of  a desired  signal  in  the  presence  of  other  signals  such 
as  coda. 

It  is  to  be  recognized,  of  course,  that  the  potential  reduction 
in  the  noise  level  by  using  the  VLP  stations  may  not  mean  that  the 
Polarization  and  Azimuthal  Filter  techniques  are  optimum  in  the  VLP 
signal  range.  The  first  reason  is  that  the  reduction  in  background 
noise  levels  for  the  VLP  instrumentation  may  be  accompanied,  at 
some  sites,  by  a reduction  in  amplitude  for  the  events  of  interest. 
Thus,  the  desired  signal  to  background  noise  level  improvement  may 
not  be  as  great  as  originally  conjectured.  The  second  reason  is  that 
the  use  of  longer  period  signals  can  result  in  more  frequency-time 
overlaps  of  the  desired  event  and  the  interference.  For  signals  having 
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periods  from  25  to  60  seconds,  the  time  aperture  for  evaluating  the 
spectral  components  is  limited  to  about  200  seconds.  This  is  the 
duration  over  which  the  multipath  propagation  is  stable10.  There  are, 
thus,  only  six  independent  frequency  bins  and  the  probability  of  the 
desired  event  and  the  interfering  event  randomly  occupying  the  same 
frequency  bin  is  1 /6th. 

When  this  simultaneous  occupancy  occurs,  the  interfering 
event,  when  large,  will  usually  suppress  the  desired  event  unless 
fortuitously  structured.  If  the  critical  duration  of  the  desired  event 

Q 

is  600  seconds,  only  ( 5/6  ) or  0.  58  or  about  half  the  time  will 
the  desired  event  remain  unobscured  by  the  interference.  If  the 
critical  duration  is  longer,  the  performance  is,  of  course,  worse. 

The  problem  is  that  the  desired  event  and  the  interference 
must  be  frequency  separable.  However,  at  least  half  the  time,  for 
VLP  signals,  this  frequency  separation  will  not  be  present  and  both 
the  Polarization  and  Azimuthal  filter  techniques  will  fail.  This  is 
the  reason  the  Ellipticity  Filter  has  been  proposed  as  a technique 
' which  will  operate  whether  or  not  the  events'  frequencies  overlap'11’  14 . 
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2.  0 Polarization  Filter 

The  Polarization  Filter,  as  developed  by  Choy  and  McCamy\ 
separates  a Rayleigh  wave  from  Love  waves  by  retaining  and 
recombining  only  those  signals  which  possess  the  desired  attribute 
of  a ninety  degree  phase  shift  between  corresponding  frequency- 
time bin  components  of  the  vertical  and  horizontal  outputs  of  v.  three- 
axis  seismometer  site.  The  specific  procedure  is  to  select  time 
windows  of  about  two-hundred  seconds  duration.  Ihis  is  long  enough 
to  have  about  a half-dozen  independent  frequency  bins  over  the  VLP 
range  of  interest  (25  to  60  seconds)  but  snort  enough  so  that  the 


multipath  is  somewhat  stable.  The  descrete  Fourier  Transforms  of 
the  selected  horizontal  and  the  vertical  seismometer  outputs  are  then 
computed  for  over-lapping,  time -shifted  windows.  Ideally,  the  selected 
horizontal  component  is  the  radial  component,  when  the  azimuth  is 
known,  or  the  corresponding  radial  component  when  a particular 
azimuth  is  being  tested.  However,  the  selection  should  not  be  too 
critical  of  azimuth  provided  the  desired  horizontal  component  is  not 
near  tangential. 

For  each  frequency  component  in  the  VLP  range,  the  phase 
difference  between  the  vertical  and  selected  horizontal  output  is 
evaluated.  Each  frequency  component  is  then  either  passed  or 
attenuated  based  on  the  value  of  this  phase  difference. 
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Rather  than  using  a pass  or  stop  criterion,  Choy  and  McCamy 
used  a filter  function  with  a monotonic  fall-off  which  depended  on 
the  absolute  value  of  the  phase  difference  between  the  vertical  and 
the  selected  horizontal  component.  The  particular  filter  function 
used  by  Choy  and  McCamy  is 

F ~ An2N  L ^CiS)  ~ ^(oS)J 

in  which  (j  (<**)  and  4^^  are  P^ases  °f  the  selected  vertical  and 

"S 

horizontal  component. 

As  will  be  seen,  the  exact  filter  function  is  relatively  unim- 
portant although  one  could  optimize  the  functional  form  depending  upon  the 
the  expected  relative  amplitudes  of  the  desired  signal,  interference, 
and  background  noise.  Since  the  basic  interest  is  in  separating  out 
a desired  signal  when  the  interference  is  high  and  the  background 
noise  level  low,  a near  rectangular  filter  is  optimum.  In  this  sense. 


tt  » 
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this  filter  function  is  equivalent  to  passing  all  frequency  components 
whose  phase  angle  differences  from  the  90°  criterion  are  less  than 
some  A and  stopping  all  frequency  components  whose  difference 
exceeds  A . For  the  filter  function  used  by  Choy  and  McCamy,  the 
pass  range  A depends  on  the  exponent  N , large  N corresponding  to 
a smaller  A . 
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The  filter  output  is  calculated  by  multiplying  each  Fourier 
component  of  the  signal  by  the  filter  function,  F . The  output 
frequency-time  bin  representations  of  the  passed  signal  are  then 
transformed  back  into  the  time  domain,  if  desired,  and  smoothed, 

The  resulting  signal  is,  thus,  composed  of  frequency-time  bin  signals 
whose  phase  differences  between  the  vertical  and  horizontal  are  near 
ninety  degrees  - as  for  a pure  Rayleigh  wave. 
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2.  1 Suggested  Improvements  in  the  Polarization  Filter 

The  suggested  improvement  in  the  Polarization  I*  ilter  is  to 
configure  an  Ellipticity  Filter.  More  modest  improvements  are 
possible  depending  on  the  specific  multipath  components  present. 


2.  1.  1 Rayleigh- Love  Wave  Mixed  Event 

The  Polarization  Filter  was  originally  configured  to 
separate  Rayleigh  from  Love  waves.  Since  Love  waves  are  horizontally 
polarized  with  no  vertical  component  of  motion,  a mixed  signal  will 
pass  only  at  those  frequencies  where  the  Rayleigh  wave  dominates, 
or  where  the  horizontal  components  of  the  Rayleigh  and  Love  signals 
are,  by  chance,  completely  in  phase  or  out  of  phase.  The  situation 

is  as  shown  below. 


Horizontal  Components 


Vertical  Components 
Rotated  90° 
for  Convenience 


What  the  filter  does  is  separate  mixed  Rayleigh  and  Love  w es  which 
do  not  occupy  the  same  frequency-time  bin  at  the  same  time. 


The  Polarization  Filter  could  be  improved  for  this  mixed 


event  situation  by: 


1,  For  each  frequency  bin,  test  the  amplitude  as  well  as  the 
phase.  Only  those  frequency  bins  in  which  there  is  both 
significant  horizontal  and  vertical  components  should  be 
considered  to  have  a Rayleigh  wave  component.  This 
would  eliminate  some  background  noise  as  well  as  the 
happenstance  in  which  the  .hove  wave  component,  albeit 
alone  in  a frequency  bin,  appeared  to  be  a Rayleigh  wave 
because  of  a near  ninety  degree  relative  phase  angle  with  a 
corresponding  vertical  noise  component. 

2.  A further  refinement  is  to  consider  only  those  frequency 
components  to  be  a Rayleigh  wave  in  which  the  ratio  of 
the  horizontal  to  vertical  component  is  within  the  allowed 
range  of  ellipticity  values  (including  measurement  errors), 
including  amplitude  as  well  as  phase. 

This  latter  refinement,  of  course,  would  bring  the  perfor- 
mance of  the  polarization  Filter  closer  to  that  of  the  Ellipticity  Filter. 
Then,  the  Polarization  Filter  would  separate  Rayleigh  from  Love  waves 
even  for  the  happenstance  situations,  provided  the  waves  did  not 
occupy  the  same  frequency-time  bin.  Complete  rejection,  and  thus 
error,  would  occur  if  any  significant  portion  of  these  waves  did  over- 
lap. The  Ellipticity  Filter  would  separate  the  signals  even  with  an 
overlap  condition. 


2.1.2  Rayleigh- Rayleigh  Wave  Mixed  Events 


While  the  Polarization  Filter  was 
Rayleigh -Rayleigh  wave  mixed  events,  it 
this  situation  to  determine  if  it  could  be 


not  designed  to  handle 
is  desirable  to  examine 
adapted  to  do  so. 
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Pure  Rayleigh  waves  have  the  characteristic  that  the  vertical 

,o 


components  are  always  i 90  out  of  phase  with  both  horizontal  com- 
ponents. If  the  two  Rayleigh  waves  do  not  overlap  in  the  frequency- 
time bins,  it  will  be  quite  evident  that  two  Rayleigh  waves  are  present. 
The  basis  of  the  sort  can  then  be: 


The  smooth  behavior  of  the  instantaneous 
frequency  with  time. 


2. 


The  smooth  behavior  of  the  instantaneous 
amplitude  with  time. 


When  two  Rayleigh  waves  are  present  in  one  frequency-time  bin, 
the  composite  signal  is  as  shown  below. 


total  vertical 


Horizontal  Components 


Vertical  Components 
Rotated  90  degrees 


V V 
i > z 


are  the  vertical  components 


are  the  ellipticities  of  the  two  Rayleigh  waves, 
since  these  may  be  bearing  dependent. 


*•> 


are  the  corresponding  bearings. 
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When  the  two  waves  are  from  the  same  azimuth,  the  total 
horizontal  component  will  be  90  degrees  with  respect  to  the  total 
vertical  component.  There  is  no  way  to  separate  this  type  of 
mixed  event  vising  a single  seismological  station.  If  the  azimuths 
differ  significantly,  however,  the  mixed  signal  will  be  identifiable 
as  not  being  a single  Rayleigh  event,  provided  the  azimuth  dif- 
ference is  not  closely  compensated  by  a change  in  ellipticity.  It 
will  not  be  possible  to  determine,  however,  if  the  mixed  event  is 
Rayleigh-Rayleigh  or  Rayleigh- Love  unless  the  time  behavior  is 
examined  and/or  the  range  of  potential  ellipticity  values  is 
considered  and/or  the  bearing  is  known. 

What  the  filter  does  do  is  to  help  separate  mixed  Rayleigh 
waves  which  do  not  occupy  the  same  frequency  bin  at  the  same  time 
and  reject  mixed  Rayleigh  waves  which  do  occupy  the  same  frequency 
bin  at  the  same  time  --  provided  they  are  not  from  the  same  azimuth. 

2.  1.  3 Love-Love  Wave  Mixed  Events 

Mixed  signals,  comprised  of  only  Love  waves,  will  not  pass 
at  all  because  the  Polarization  Filter  is  dependent  on  the  phase 
difference  between  vertical  and  horizontal  components.  Since  Love 
waves  are  horizontally  polarized  and  have  no  vertical  component, 
the  filter  will  reject  the  entire  record  unless  the  vertical  back- 
ground noise  fortuitously  resembles  the  equivalent  vertical  Rayleigh 
component.  A direct  amplitude  ter,  should  resolve  this  situation. 
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2.  2 Performance  Evaluation  of  the  Polarization  Filter 

The  Polarization  Filter  is  a technique  for  sorting  out  Rayleigh 
waves  in  a mixed  event  situation.  The  basis  of  the  sort  is  that  the 
horizontal  and  vertical  components  of  a Rayleigh  wave  have  a ninety 
degree  phase  difference.  When  interference  is  present,  the  corres- 
ponding phase  difference  is  usually  not  ninety  degrees  so  that  the 
presence  of  the  interference  can  be  recognized. 

In  the  Polarization  Filter,  the  two  input  signals  are  Fourier 
transformed,  using  a high  speed  algorithm  such  as  the  FFT,  and 
the  complex  outputs  of  the  corresponding  frequency  and  phase  compar 
In  the  absence  of  interference,  the  desired  signals  are  perturbed  by 
the  additive  background  noise  to  give  the  total  observed  signals.  Thu 
the  two  observed  signals  are  related  to  the  desired  signal  by 


mTO»*«  ' H^-vra^S^OW  wo*  *»f 


The  numerator  is 


^e  R|<,  ^<j  **^  N^fjio)  - Jh i R^  (^ oS)  ge  N^^to) 


so  that 


A <PK 


cos  <PK  5m  ^kQw)  ~ sin  ^ £e  ^K  ^w) 

*K 


in  which  is  the  amplitude  of  the  desired  component. 

One  notes  that,  for  Gaussian  noise,  the  inphase  and  quadrature 
noise  components,  £eA!K^uj)  , jv*  fslK  ^ uj)  are  Gaussian.  Hence, 
the  distribution  of  the  phase  angle  error  is  Gaussian  since  it  is  the 
sum  or  difference  of  two  weighted  Gaussian  distributions. 

The  mean  value  of  the  Gaussian  distribution  is  zero,  since 
the  expected  value  (denoted  by  brackets)  is 


-s  A > 


cos  h<i  bn  Nj^  Quj  )>  “ sin  efK  -v  Re  N ^ (it**)  ^ 


The  variance  is 


V - <(*<pKr>  = 


cos2yK<[3rr,NK(j,'J)J''/  ).  sin^^Jf  ^[Ve  N|<6 


K 


K 


since  the  crossterm  is  zero  - i.  e.  , ^ A ^ ~ O. 
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If  the  noise  has  a level,  <5^ 


^[J„NK(^)JS>  = <l^HK(r)Y>  = <r* 


k ‘<r 


hK 


then 


/r-  2 


= ri 


- 2. 


hK 


/ * 
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Thus,  the  error  in  one  phase  angle  is  Gaussian  with  zero 


mean  and  variance 


crft; 
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R, 


= n sv' 


For  the  other  phase  angle,  the  error  is  also  Gaussian  with  zero 
mean,  but  variance 


‘ H2 


H 


K 


nsG 


The  variances  of  the  phase  angle  error  are  thus  proportional  to  the 
noise -to -signal  ratios  (nsr)  when  the  noise  levels  are  small.  The 
difference  angle,  which  is  the  argument  of  the  filter  function,  is 
thus  normal  with  zero  mean  as  well.  When  the  snr  in  both  seis- 
mometer signals  are  about  the  same,  the  variance  is  twice  the  nsr 
of  one  of  the  inputs;  when  the  snr  is  much  smaller  for  one  seis- 
mometer signal  than  is  the  other,  the  variance  is  the  large,!  noise- 

to-signal  ratio.  In  general 


<? 


- h SKj  +-  ns 


for  the  difference  angle. 
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The  filter  function,  F(jiO}(f 0 , can  have  a variety  of  forms, 

one  being  that  selected  by  Choy  and  McCamy. 

F , <p")  - sin™  ( <Pf  - (py) 

However,  the  basic  effect  of  the  filter  function  is  to  pass  all  sig- 
nals for  which 

/ 0,  - Qx  I - a 

and  inhibit  all  other  signals.  The  performance  of  the  system  can 
thus  be  adequately  described  by: 


1,  An  equivalent  rectangular  filter  function,  + A wide 

2.  The  operating  characteristics  in  terms  of  the  pro- 
bability of  accepting  a desired  signal  p(A/D)  vs 
the  probability  of  accepting  interference  p(A/I). 


This  characterization  will  be  done  in  the  next  two  sections. 


1 

k « 


n 


I ! ! 


2.2.1  The  Equivalent  Rectangular  Filter  Function 


There  are  several  ways  an  equivalent  rectangular  filter 
function  can  be  defined.  The  most  rigorous  definitions  depend  on 
the  statistics  of  the  interference  to  desired  signal  level  and  the 
attenuation  required  to  render  the  interference  significantly  below 
the  desired  signal  level.  This  approach  is  a little  too  involved 
for  the  general  considerations  of  interest  here. 

Instead,  three  reasonable  definitions  will  be  considered  and 
the  results  shown  to  be  essentially  the  same. 


1.  RMS  Window 

2.  Integrated  Shape  Window 

3.  Six  dB  Level  Window 


The  equivalent  rectangular  filter  function  can  be  based  on  an 


1 rms  window"  defined  as 


-rC. 


J & 2 6 d& 


t r/x 


/ cos2*  6 d& 


Since  lfjz 
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cos  d<9 
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the  equivalent  A can  be  solved  for  by  iteration. 


With  equal  intuition,  an  integrated  shape  window  can  be  defined 


J"  Cos^  & d& 


so  that  the  areas  are  equal,  Then 


i.  3.  5"...  (2H  - l)  V 


2.4.fc...  (*H)  2 


With  equal  intuition,  a third  definition  can  be  used,  the 


6 dB  level  window 


Cos 2 A = >/Z 


-'('/>) 


The  results  are  shown  in  Table  I.  The  different  definitions  do 


not  lead  to  significantly  different  results. 


Experimentally,  N = 4 was  found  satisfactory  fo**  the  few  evaluations 


made  with  data  . This  corresponds  to  passing  signals  whose  phase 


difference  between  the  vertical  and  horizontal  is  within  +20  degrees 


of  90  degrees.  One  also  notes  that  the  effective  aperture  is  not  a 


strong  function  of  N for  large  h/  . 
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Table  I:  Equivalent  Rectangular  Filters 

for  the  Power  Law  Sine  Window 


N 

RMS 

Integrated 

6 dB 

A radians 

A degrees 

Shape  A 

A 

1 

0. 5678 

32.  53 

45.  00° 

45.  00° 

2 

0. 4444 

25.46 

33.  75° 

32.  77° 

3 

0. 3767 

21.  58 

28. 13° 

0 

27.  01 

4 

0.  3327 

19.  06 

24.61° 

23.  51° 

5 

0.  3011 

17.  25 

22. 15° 

21.  09° 

A is  a ’ half"  window  definition 
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2.  2.  2 Operating  Characteristics 


The  probability  of  accepting  the  desired  signal  is 


+ A -<Pyo(hst'l  -f  nsr  ) 
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Wfl/lO  = f , i 

-A  A 2T\  [ t\srl  -t-  hsrA 


n'/z 


- X' 
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A/Yhsr,  fhsi^)  ^ 


J 2 TTn 


cl  X 
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= 2 
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p 
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A/fnsr;  +hff«'a 


y/A  a/^F 


= I - 2 erfc 


^ ^hsr(  +•  hs r^)  ^ 1 


in  which  £rfc  is  the  complement  of  the  error  function 
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The  interference,  when  present,  is  assumed  to  be  much  larger 
than  the  signal.  The  corresponding  phase  angle  <P  then  has  almost 
a uniform  distribution.  The  probability  of  improperly  accepting  the 
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interference  is  then 


?(/)/£)  = zj 


These  two  expressions  give  the  operating  characteristics  of  the 


Polarization  Filter 


T (e\  /rO 


f 1 “ 2 erfc  ( (nsr<  Asrp»>  ^ 


PO/l) 


■ ['  - 2 e*fc  [prp o'/*’]] 


The  difficulty  in  obtaining  satisfactory  performance  from  the 


Polarization  Filter  can  be  seen  from  Figure  1.  Signal-to-noise 


ratios  of  about  20  dB  are  required  for  satisfactory  performance. 


The  key  advantage  of  the  Polarization  Filter  can  also  be  seen.  The 


non-linear  processing  will  reject  very  high  levels  of  interference. 


This  performance  threshold  can  be  described  analytically. 


When  T(f\/X)  is  to  be  small 

„ 7 r P(A  /r) 


tP  ?C*l  f) 


^hSK  + nsO^ 
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Probability  of  Accepting  the  Desired  Event,  P(A/D) 
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in  which 
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i nsri 

= noise-to- signal 

ratio 
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for  vertical  component  in 

a common  frequency  bin 

nsr 

•••  noise-to- signal 

ratio 

for  horizontal  component 

in  a common  frequency  bin 
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and 


P(a  ] p)  = JIt? 


?(*/!) 


(nSKt  +■  hsrj) 
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Thus,  to  simultaneously  have  a high  acceptance  of  the  desired  signal 
and  a low  acceptance  of  the  interference  requires 


?^/3>) 
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and 


?(n/i) 


( t\Srt  +•  hs^) 
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For  a reasonable  level  of  accepting  interference  of  ten  percent, 

^nsr(  +■  ftsra)  = 0.  10  J ZTT 

= 0.251 


or 


(hsrt  + hSK^)  = 0.063 

10 log  ^hsr(  f hsra)  = -12  dB 

Thus,  the  ratio  of  the  desired  signal  to  the  background  noise  (in  the 
frequency  bin  which  has  the  desired  signal)  should  be  12  to  15  dB  or 
better  for  good  performance.  Since  there  are  about  a half-dozen 
frequency  bands  over  the  frequency  range  of  interest,  the  wideband 
snr  must  be  about  4 to  7 dB  for  good  performance.  If  the  reason- 
able level  of  accepting  interference  were  set  at  five  percent 
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The 
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hs k(  + hS  Kj')  = o.  016 


10  log  6isk  + hsr  ~)  = -18  dB 


The  requirement  on  the  ratio  of  the  desired  signal  to  background 
noise  level  would  then  be  18  to  21  dB.  The  corresponding  wide- 
band snr  should  then  be  10  to  13  dB  for  good  performance. 

The  remaining  point  is  to  note  that,  for  the  level  of  accept- 
ing interference  of  ten  percent,  the  window  size  is  given  by 


? I a / T) 


i/r 


0.  10 


or  the  window  is  + 18°.  Thus,  a ten  percent  interference  accept- 
ance level  corresponds  to  the  phase  window  of  ± 20°  which  was  used 
by  Choy  and  McCamy.  Note  also  that  a of  10  percent. 


on 


a per  frequency  bin  basis,  means  for  six  bins  no  interference  will 

\6 


be  passed  with  a probability  of  (0.9)  or  only  fifty  percent. 

This  means  that,  with  a desired  signal  to  background  noise 
level  of  4 to  7 dB,  the  Po  .<  -ization  Filter  would  operate  satisfactorily 
against  high  interference  levels  about  half  the  time.  The  filter 
function,  however,  should  be  more  rectangular  than  that  used  by 
Choy  and  McCamy  and  with  a window  of  t 20°.  To  improve  the 
level  of  satisfactory  performance  to  about  seventy-five  percent  of 
the  time,  the  desired  signal  to  background  noise  level  would  have  to 
be  10  to  13  dB,  and  the  filter  function  more  rectangular  and  with  a 
window  of  ± 10  degrees. 
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3.  0 Azimuthal  Filter 

The  Azimuthal  Filter  of  Choy  and  McCamy  is  designed  to 
separate  Love  waves  from  interfering  waves.  Love  waves  are 
horizontally  polarized  transverse  waves  i.  e.  , the  particle  motion 
is  perpendicular  to  the  direction  of  propagation  and  exists  only  in 
the  horizontal  plane.  The  basis  of  the  sorting  procedure  is  to  deter- 
mine the  apparent  azimuth  direction  from  the  horizontal  components 
and  to  pass  only  those  waves  whose  apparent  azimuth  direction  is 
close  to  the  azimuth  of  the  desired  event.  Filters  of  this  type  have 
been  described  by  Sax  and  Mims  (19&5)  and  Simons  (19&8).  Corres- 
ponding counterparts  are  used  in  radio  direction  finding  and  underwate 
sound  systems. 

The  processing  system  is  very  similar  to  that  of  the  Polari- 
zation Filter.  The  N - S component  and  the  E - W component  of 
the  seismometers  are  Fourier  transformed.  For  each  frequency 
bin,  an  apparent  azimuth  is  calculated 


X£w 

X„5  £j-> 


is  the  E - W component  in  the  frequency  bin  u» 
is  the  N - S component  in  the  frequency  bin  uj  . 


A 

The  difference  between  the  apparent  azimuth,  B , and  the  azimuth 
of  the  desired  event,  B , (or  the  azimuth  being  tested) 
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determines  whether  or  not  the  frequency  components  are  passed. 

When 

0M  = ± ^o°  t A 

the  Love  wave  is  passed.  Otherwise  the  frequency  component  is 
attenuated  or  blocked. 

Choy  and  McCamy  actually  used  the  filter  function 

F[J)  = Sin  ^ 6 (uj) 

as  the  blocking  function.  As  analyzed  for  the  Polarization  Filter, 
the  exact  filter  function  is  not  too  critical,  although  one  could 
optimize  the  function  form  for  particular  levels  of  desired  signal, 
interference,  and  noise.  The  filter  function  essentially  operates 
by  passing  all  signals  whose  azimuths  are  within  some  i A of  the 
expected  azimuth.  The  relationship  between  the  exponent,  N,  of  the 
cosine  power  law  filter  and  A is  the  same  as  that  given  in 
Section  2.  2.  1. 

As  with  the  Polarization  Filter,  the  filter  output  is  calculated 
by  multiplying  each  Fourier  component  of  the  signal  by  the  filter 
function.  The  frequency-time  bin  representations  of  the  passed  signal 
are  then  transformed  into  the  time  domain  and  smoothed  to  give  the 
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estimated  Love  wave  output  for  signals  coming  from  the  desired 
azimuth. 

The  limitations  on  this  filter  are  similar  to  those  of  the 


Polarization  Filter. 


A high  signal-to-background  noise  level  is 
required. 


The  filter  sorts  rather  than  filters,  so  that  the 
the  desired  signal  and  the  interference  should 
not  occupy  the  same  frequency  bins  at  the  same 
time. 
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3*  1 Suggested  Improvements  in  the  Azimuthal  Filter 


The  Ellipticity  Filter  is  the  suggested  improvement  of  the 
Azimuthal  Filter  since  it  combines  both  the  features  of  the  Polari- 


zation Filter  and  the  Azimuthal  Filter  as  well  as  other  decided 
advantages.  More  modest  improvements  can  be  made  depending  on 
the  specific  multipath  components  present. 


3.  1.  1 Love-Love  Wave  Mixed  Events 


The  Azimuthal  Filter,  as  configured,  works  as  expected 
against  Love-Love  wave  mixed  events.  The  system  will  filter  Love 
waves  from  Love  waves  having  a different  azimuth.  Confusion  will 
result  only  if  the  mixed  event  is  due  to  Love  waves  coming  from  the 
same  azimuth  or  from  directions  180°  apart.  A simple  improvement 
is  to  test  amplitude  as  well  as  azimuth  so  that  background  noise  levels 
in  those  frequency-time  bins  which  are  void  of  either  Love  wave 
components  are  not  confused  with  the  Love  wave  components. 


3.1.2  Love -Rayleigh  Wave  Mixed  Events 


For  Love-Rayleigh  wave  mixed  events,  the  Love  wave  will 
be  properly  selected  unless  the  Rayleigh  wave  occupies  the  same 
frequency-time  bins  as  the  Love  wave  or  the  Rayleigh  wave  coming 
from  an  azimuth  which  is  normal  to  the  expected  Love  wave  azimuth. 
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In  this  latter  case,  examination  of  the  frequency  of  the  vertical  com- 
ponent of  the  Rayleigh  wave  will  serve  to  indicate  which  is  the  Love 
wave  and  which  is  the  Rayleigh  wave. 


3.  1.3  Rayleigh -Rayleigh  Wave  Mixed  Events 


If  the  apparent  direction  of  the  composite  Rayleigh-Rayleigh 
wave  mixed  event  is  normal  to  the  expected  Love  wave  direction, 
the  system  will  confuse  the  two.  The  obvious  solution  is  to  use  the 
vertical  component  to  determine  the  presence  of  two  Rayleigh  waves. 
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3.  2 Performance  Evaluation  of  the  Azimuthal  Filter 


The  performance  evaluation  of  the  Azimuthal  Filter  is  similar 
to  that  of  the  Polarization  Filter.  The  main  difference  is  that  the 
evaluation  of  azimuth  cannot  be  done  to  the  proper  quadrant,  where- 
as the  evaluation  of  the  polarization  angle  can  be  done  to  the  proper 
quadrant.  The  result  is  that  an  acceptance  angle  of  ± A for  the 
Azimuthal  Filter  means  that  interference  will  be  effective  over  four 
regions  as  shown  in  Figure  2.  For  a reasonable  A of  20° 

(see  Section  2.  2 of  the  Polarization  Filter  analysis)  this  means  that 

there  is  an  80°  range  over  which  Love  wave  interference  is  serious 
o 

and  an  80  range  over  which  Rayleigh  wave  interference  can  be 
serious  when  the  interference  is  in  the  same  f requency-time  bin 
as  the  desired  signal.  Clearly,  the  advantage  of  the  Azimuthal 
Filter  is  similar  to  that  of  the  Polarization  Filter  --  namely,  that 
the  filter  shows  when  the  desired  signal  alone  occupies  a frequency- 
time bin  and  when  the  desired  signal  is  contaminated  by  an  inter- 
ference in  the  same  frequency-time  bin.  Neither  filters  against 
this  interference  condition,  they  merely  indicate  its  presence.  The 
Ellipticity  Filter,  on  the  other  hand,  actually  filters  against  these 
frequency-time  bin  overlapping  events. 
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SECTION  II 
MEASUREMENT  OF 

SEISMOMETER  MAGNIFICATION  AND  ELLIPTICITY 
AT  A VLP  STATION 


As  a possible  solution  to  the  mixed  event  problem  in 
nuclear  seismic  discrimination,  OAS  has  developed 
a filter  termed  the  "ellipticity"  filter  which  has  the 
superposed  Rayleigh  and  Love  waves.  Briefly,  the 
ellipticity  filter  processes  the  output  signals  of  a three- 
component  seismometer  station,  using  the  value  of 
ellipticity  at  the  station  to  produce  the  complex  am- 
plitudes and  azimuths  of  the  superposed  surface  waves. 

In  order  to  use  this  filter,  both  the  station  calibration 
and  ellipticity  must  be  known.  In  this  report,  the 
calibration  was  determined  for  the  Kongsberg  VLP 
station  during  March  1972.  With  this  calibration  and 
three  (3)  earthquake  seismograms,  the  value  of  sta- 
tion ellipticity  was  estimated.  The  rms  error  is  less 
than  6%  for  periods  from  21  to  64  seconds. 


1.  0 Introduction 

As  a possible  solution  to  the  mixed  event  problem  in  nuclear 

seismic  discrimination,  OAS  has  developed  a filter  termed  the 

"ellipticity”  filter  which  has  the  capability  of  separating  sets  of 

11  12 

superposed  Rayleigh  and  Love  waves  ' . Briefly,  the  ellipticity 

filter  processes  the  output  signals  of  a three  component  seismometer 
station  using  the  value  of  ellipticity  at  the  station  to  produce  the 
complex  amplitudes  and  azimuths  of  the  superposed  surface  waves. 

In  order  to  use  this  filter,  both  the  station  calibration  and 
ellipticity  must  be  known.  In  this  report,  the  calibration  is  deter- 
mined for  the  Kongsberg  Y.LP  station  during  March  1972.  With 
this  calibration  a:.d  the  earthquake  seismograms,  the  value  of 
station  ellipticity  is  estimated.  The  mis  error  is  less  than  6% 
for  periods  from  2 to  64  seconds. 

Another  approach  to  the  determination  of  ellipticity  could 
have  been  a theoretical  calculation  based  on  the  known  crustal 
structure  at  Kongsberg,  Norway14.  However,  this  approach  was 
not  implemented  because  the  ellipticity  measurement  program 
allowed  the  evaluation  of  the  concept  of  ellipticity  on  a group  of 
events  which  were  not  highly  selected,  and  in  addition,  permitted 
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a verification  of  the  calibration  method.  Results  of  the  calibration 
and  ellipticity  measurement  indicate  that,  with  slightly  better  instru- 
mentation than  available  at  the  time,  the  actual  separation  of  mixed 
events  could  be  achieved. 

The  calibration  was  a simple,  but  tedious  procedure.  The 
Kongsberg  station  was  selected  because  of  its  proximity  to  nuclear 
testing  sites  and  because  of  the  availability  of  digital  data  tapes. 

Ti  e rectangular  pulses  supplied  daily  to  the  system  were  used  for 
calibration.  Calibration  was  done  by  a complex  amplitude  fit  in 
the  frequency  domain.  Since  it  could  not  be  determined  whether 
the  variation  in  the  spectrum  of  the  pulses  was  due  to  change  in 
the  seismometer  parameters  or  to  noise,  four  different  calibration 
procedures  were  evaluated: 


1.  Fitting  the  daily  response  to  a triple  tuned  model, 

2.  Fitting  the  average  daily  response  to  a triple  tuned 
model, 

3.  Fitting  the  daily  response  by  parabolic  curve  fitting, 

4.  Fitting  the  average  daily  response  by  parabolic  curve 
fitting. 


The  evaluation  involved  the  use  of  each  calibration  to  determine  sta- 
tion ellipticity  from  six  events.  The  criterion  for  acceptance  of  a 
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procedure  was  the  most  consistent  ellipticity  evaluated.  In  the 
actual  estimation  of  ellipticity,  only  events  which  had  azimuths 
close  to  a cardinal  direction  were  used.  For  these  events,  cali- 
bration errors  cause  only  a scaling  error  in  the  estimation  and 
not  an  error  in  the  Rayleigh-Love  wave  separation  process. 

By  the  fitting  of  a 15  day  average  of  the  Fourier  transforms 
of  pulse  to  a triple-tuned  seismometer  response  model,  the  best 
values  (Table  II)  for  the  magnitude  of  the  complex  ellipticity  were 
obtained.  However,  the  best  values  for  the  angle  of  ellipticity 
were  obtained  by  parabolic  curve  fitting  of  both  the  amplitude  and 
phase  to  the  15  day  average  of  the  Fourier  transforms.  The 
results  are  shown  in  Table  III. 

Section  Z contains  a brief  review  of  the  VLP  System.  In 


!l 

Section  3,  a detailed 

♦ 4 

appears.  Section  4 

ii 
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for  six  earthquakes. 
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System  De scrip;  ion 

The  high  gain,  long  period  station  developed  by  Lamont- 
Doherty  Geological  Observatory  consists  of  three  seismometers, 
one  vertical  and  two  ho rizontal^  mounted  at  right  angles  to  each 
other.  Figure  3 shows  a block  diagram  of  the  station.  . Each 
seismometer  has  two  velocity  transducers,  a displacement  trans- 
ducer, a primary  calibration  coil,  and  a secondary  calibration 
coil.  In  normal  station  usage,  only  the  primary  coil  is  used  for 
calibration. 

For  the  high  gain  velocity  channel,  a velocity  transducer  is 
connected  to  a galvanometer  with  an  attached  mirror.  After  the 
position  of  the  light  beam  reflected  by  this  mirror  is  sensed  with 
a pair  of  phototubes,  the  resulting  signal  is  amplified  and  filteiod. 
After  the  signal  is  digitized  at  a rate  of  512  samples  per  second, 
it  is  averaged  for  one  second  and  then  stored  on  magnetic  tape.  In 
addition,  it  is  recorded  photographically  on  a drum  recorder.  For 
the  low  gain  velocity  channel,  the  output  of  the  other  velocity  trans- 
ducer is  recorded  directly  on  a photographic  drum  with  a recording 
galvanometer. 

In  the  displacement  channel,  a variable  capacitance  transducer 
together  with  an  oscillator  and  discriminator  are  used  to  sense  the 
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position  of  the  seismometer  boom.  The  signal  is  then  low  pass 
filtered,  digitized,  and  averaged.  The  displacement  data  is  recorded 
on  magnetic  tape  once  every  five  seconds.  In  addition  to  giving  an 
indication  of  boom  position,  the  displacement  channel  monitors 
earth  motion  from  a period  of  about  10  seconds  to  D.  C. 
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3.  0 Calibration 


In  order  to  calculate  the  station  ellipticity  from  the  digitized 
high  gain  velocity  outputs  of  the  three  component  seismometer,  the 
instrumental  response  of  the  seismometer  must  first  be  removed 
from  the  signals.  Let  Xltt,  and  ?(£)  be  the  E/\Y,  N/S,  and 

vertical  component  signals  expressed  in  digital  counts.  Assume  that 
a Fourier  transform  of  these  signals  has  been  taken  for  some  desired 
window  and  are  expressed  as  Y;U.  Y;(  jw)  , and  2: ; i'yJ)  . Then 
the  Fourier  transform  of  the  actual  ground  motions 
and  may  be  expressed  as 


y.-  (jw) 


Y;  (y*) 


Zatjw') 


X;  (jA 


where  tAytjw)  . and  w'l  are  the  transfer  functions 

of  the  component  seismometers  expressed  in  counts /micron.  The 


transfer  function  is  also  called  the  magnification  of  the  seismometer. 
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Since  the  Fourier  transforms  of  the  digitized  outputs  of  the  seismometers 
are  expressed  in  counts,  the  transforms  of  the  calculated  ground  motions 
will  be  microns. 


The  transfer  function  of  a component  seismometer  is  determined 
by  passing  a known  current  through  its  primary  calibration  coil  and 
dividing  the  Fourier  transform  of  its  output  by  the  transform  of  the 
known  input  signal.  Corresponding  to  a force  applied  to  the  seis- 
mometer, the  current  causes  a relative  acceleration  of  the  seis- 
mometer in  the  same  manner  that  earth  motion  does.  Before  this 
method  can  be  used  at  a station,  the  electromechanical  constant,  G, 
for  the  coil  must  first  be  determined.  A variation  of  the  standard 
weight  lift  method  is  used.  In  tins  method,  a weight  is  attached  to 
the  seismometer.  This  weight  causes  a static  deflection  of  the  seis- 
mometer. A D.  C.  current  is  then  sent  through  the  coil  and  its  mag- 
nitude is  adjusted  so  that  the  seismometer  returns  to  its  initial  position. 
The  value  of  G may  then  be  determined  from  the  formula 


G = 


applied  force  in  newtons 


restoring  current  in  amperes 


One  known  signal  input  for  calibration  is  a set  of  sine 
waves  with  periods  of  30,  40,  70,  and  100  seconds.  However, 
for  the  month  selected,  no  sine  wave  calibration  was  available 
in  Kongsberg.  In  1972,  the  time  between  sin-  wave  calibrations 
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at  tlie  station  was  greater  than  six  (6)  months.  Inspection  of 
Alberquerque  Seismological  Center  records  showed  calibration 


changes  of  about  10%  a month.  Thus,  sine  wave  calibrations 
could  not  be  used  even  if  they  were  present. 

In  addition  to  the  sine  waves,  rectangular  pulses  of  cxirrent 
are  applied  to  calibration  coils.  At  Kongsberg,  calibrations  pulses 
of  4>ua  and  50,U.a  are  applied  simultaneously  each  day  to  all  three 
seismometers.  Examination  of  the  plotted  digital  data  shows  that 
the  pulses  last  for  240  seconds.  Ideally,  the  large  50u.a  calibra- 
tion pxilsc  should  be  used  to  minimize  the  effect  of  noise  on  the 
calibration.  However,  the  output  shows  an  asymmetry  because  the 


large  pulse  makes  the  high  gain  velocity  channel  act  in  its  non- 
linear region.  Figure  4 shows  an  asymmetry  of  4 3 jo  into  the  y 
channel.  Thus,  the  small  pulses  were  used.  In  order  to  get 
frequency  domain  information  from  the  pulse,  a Dl*  1 transform 
with  a window  length  of  480  seconds  is  used.  This  length  window, 
b addition  to  letting  the  transient  response  of  the  system  to  go  to 
zero,  concentrates  all  of  the  signal  energy  in  the  odd  harmonies  of 
the  DFT  and  gives  maximum  noise  immunity. 


3.  1 Derivation  of  Calibration  Equation 


To  use  the  pulses  described  in  Section  3.  0 as  a calibration 
signal,  the  relation  between  the  amplitude  of  the  pulse  and  the 
transfer  function  must  be  derived.  Since  a current  corresponds  to 


a force  applied  to  the  seismometer 


Y 


G; 


me 


whe  re 


Y = 
G = 


seismometer  frame  acceleration  due  to 
earth  motion 


motor  constant  of  calibration  coil 


l = 
»vve  - 


current  in  calibration  coil 


effective  mass  of  seismometer, 


Because  the  force  generated  by  the  coil  is  not  applied  at  the  center 
of  mass  of  the  seismometer,  the  effective  seismometer  mass  is  re- 
lated to  the  true  mass,  Ivl,  by 


VY\0  — 


TA  icw 


whe  re 


Cm 


distance  from  the  hinge  of  the  pendulum  to 
its  center  of  mass 


rc  - distance  from  hinge  to  center  of  coil 


Therefore, 


Y 


G i rc 


M Y. 


C*r\ 
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Taking  the  Fourier  transform  of  the  equation 


n(4'a))  = 


whe  re 


V ( jw)  is  the  Fourier  transform  of  the  displacement 
amplitude  of  earth  vibration 

*•  the  Fourier  transform  of  the  input  current. 

The  transfer  function  or  dynamic  magnification,  1M  C\J) , of  the  seis- 
mometer system  is  then  given  by 


MtjvJ , 

= 


Fourier  transform  of  Digital  Record  Amplitude 
Fourier  transform  of  the  displacement  amplitude 


V r 

Y (.  \w) 


Thus 


K (\ ^ - _ Yr  Wlvfr  H rc 

All  that  remains  to  be  calculated  is  X Assume  that 


Ut)  = X cal  0 <t  ^ T 

O t>o 


Then  the  Fourier  transform  of  the  input  is 


Tt^uA  = ^ 1(F)  e-*1^  dfc 


-oD 


48 


Integration  gives 


Xtj«>  = 


2.Xcai 


V/ 


J$kf< u v.  ^ 


t'4?) 


Hence, 


2 G-Ic^  /Oaw  0=£Vc 
If  a DFT  of  length  2T  is  used,  then  the  equivalent  formula  is 


v^X 

3U 


H(jvon1=  _ Yr(hAp  w^3  Hirc>w£ 


who  re 


■Ha 

**r 


Simplifying  yields 


j Yr  HrCt 

G X ca-\  y^. 


for  n odd 


The  expression  is  not  defined  for  even  harmonics  because  the 
amplitude  of  the  input  at  even  harmonic  is  Kero. 

The  VLP  seismometer  has  the  following  constants 


Ye  rtical 

Horizontal 

13. 2 kg 

10.  0 kg 

Y^fy, 

. 230  meters 

.227  meters 

rc 

.378  meters 

.378  mete  rs 
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Using  these  constants,  the  magnifications  , VAxCjvJ*'). 

and  of  the  k,  x,  and  y channels,  respectively,  may  be 

reduced  for  odd  to 


k.oi'S  *tYrx( S^^r.3 
'G*X  ca\ 

. 3.00 3..  j 

GXclf 

3.003  jYrai(iw^'*Jn3_ 

J&xa\ 

where  and 

Yy*2  are  the  DFT  transforms  of  the  x,  y 

and  h seismometer  outputs  (the  length  of  the  transform  is  twice 
the  length  of  the  pulse). 

3.  2 Averaging  and  Curve  Fitting  of  Calibration 

In  order  to  obtain  reliable  estimates  of  the  seismometer 
magnification  from  the  small  calibration  pulses,  aveiaging  and 
frequency  smoothing  were  applied  to  the  harmonic  data  obtained 
from  the  magnification  formulas  in  the  preceding  section.  Figure  5 
shows  the  actual  measured  values  of  magnification  amplitude  and 
fitted  curves  for  three  consecutive  days.  By  just  examining  the 
curves,  it  is  impossible  to  tell  whether  the  12%  variation  in  the 
peaks  of  the  curves  should  be  attributed  to  an  actual  variation  in 
the  transfer  function  or  just  to  a random  variation  caused  by  noise. 
Therefore,  for  the  calculation  of  cllipticity,  it  was  decided  to  do 
the  calibration  on  both  a daily  basis  and  on  an  average  of  15 


days,  so  that  variations  of  the  instrument  response,  if  they  occurred, 
would  not  be  averaged  out.  The  best  method  would  then  be  the  one 
which  resulted  in  the  most  consistent  ellipticity  calculation.  As 
explained  in  the  following  paragraphs,  two  different  curve  fitting 
methods  were,  used  on  both  the  daily  and  average  calibrations. 

One  method  of  smoothing  was  to  fit  a triple  tuned  model 
to  the  amplitude  of  the  magnification.  From  the  parameters  deter- 
mined by  the  curve  fit,  the  phase  response  was  calculated.  This 
method  is  explained  in  detail  in  Appendix  A.  Figures  6 and  7 
show  the  application  of  this  method  to  both  daily  and  average 
calibrations.  Since  the  phase  response  is  determined  from  the 
amplitude  response,  the  pulses  do  not  have  to  be  precisely  aligned 
during  the  averaging  process.  Tlie  main  disadvantage  of  this 
method  is  that  even  though  it  may  give  a reasonable  fit  to  amplitude 
response,  there  may  be  systematic  error’s  in  the  phase  response. 

If  the  actual  parametric  form  of  the  transfer  function  is  not  known 
and  must  be  deduced  from  the  data,  then  the  amplitude  response 
must  be  known  very  accurately  in  order*  to  get  the  phase  response 
from  it. 

In  order  to  avoid  the  possible  phase  response  problems  and 
distortions  introduced  by  the  tuned  circuit  method,  a second  set  of 


Seismometer  Amplitude  Response  Normalized 
to  Unity  from  Triple  Tuned  Model,  Kongs- 
berg  Z Channel,  Average  Days  67  to  81 

Damping  =1.6 

Center  Period  = 45  sec. 


Period  (seconds) 


(Normalized  Amplitude  Response  ( 1.  0 corresponds  to  1145  counts /micron) 
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Fig.  7:  Seismometer  Amplitude  Response  Normalized 

to  Unity  Kongsberg  Z Channel,  Day  75 

Triple  Tuned  Model,  Damping  = 1.70, 

Center  Period  = 35  sec. 

Parabolic  Curve  Fitting 
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Normalized  Amplitude  Response  (1.0  corresponds  to  858  counts/micron) 
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both  daily  and  average  calibrations  were  derived  by  just  the  free 


hand,  parabolic  curve  fitting  of  both  amplitude  and  phase  of  the 
seismometer  magnification.  Figures  7 and  8 show  typical  results 
of  this  method. 


i 

1 

I 

; 

* 

! 

i 


...  . 


55 


r^i««vMi!  <*rW»'«*W*  -« 


4.  0 Measurement  oi  Ellipticity 


l 

I 

I 


By  the  use  of  a filter  similar  to  the  Rayleigh- Love  ellipticity 
filter  described  in  the  First  Quarterly  Repor*t  (OAS  TR  74-232)**, 
ellipticity  measurements  were  made  on  six  (6)  earthquakes  with 
four  (4)  different  calibrations.  Table  IV  contains  a list  of  events. 
Three  of  the  events  were  from  the  Honshu  region  and  three  from 
the  Aleutians.  Since  the  asstunptions  of  this  ellipticity  filter  are 
slightly  different  than  the  earlier  one,  a detailed  derivation  will 
now  be  given. 

4.  1 Description  of  Ellipticity  Measurement  Filter 

Assume  we  have  a three -component  seismometer  with 
the  output  signals  X(0  . Y(0,  and  2.60  and  that  the  Fourier  trans- 
form of  these  signals  has  been  taken.  The  signals  X,  Y.  and  5L 
represent  one  of  the  Fourier  components  of  the  transformed  signal 
and  and  represent  the  complex  magnification  or  trans- 

fer function  at  the  frequency  of  the  Fourier  components  of  the  signals. 

Under  these  assumptions,  a Rayleigh  and  Love  wave  coming 
from  a single  azimuth  may  be  represented  as 


X/Mx  = - ER  sin  6 - U cos  © (£) 

COS  0 4-  L s*n  Q (z) 


x/*\ 


^3) 


is  the  complex  amplitude  of  the  Rayleigh  wave 

is  the  complex  amplitude  of  the  Love  wave 

is  the  azimuth  of  both  the  Rayleigh  and  Love  wave 

is  the  complex  station  ellipticity  (ratio  of  horizontal 

to  vertical) 

Assuming  that  the  Love  and  Rayleigh  wave  arrive  by  the  great 
circle  path  (Second  Quarterly  Report,  OAS  TR  74-252),  the  azimuth  of 
the  event  may  be  calculated  from  data  supplied  in  the  PDE.  The  above 
three  complex  equations  represent  six  nonlinear  equations  in  the  unknown 
real  and  imaginary  parts  of  H , R , and  \ — . Thus  we  have  six 

equations  in  six  unknowns. 

Substituting  (3)  into  (1)  and  (2)  gives 

(£;)  = -E(jt)£,n9  ' u) 

( ) ° ~ ^ ( y;  ) cc,s  6 + L 501 e 

J z 


Multiply  (4)  by  Sm  0 and  (5)  by  co$(9  , and  then  adding  the  resulting 
equations  gives 


whe  re 

R 

\_ 

e 

E 
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In  a similar  manner,  the  solution  for  the  complex  amplitude 


of  the  Love  wave  is 


x-u  vv-  nuve  to  / \/  s / \/  \ 

~ (~n7)c°sB  + (~i^)sine 

(-t) 

The  value  of  the  complex  amplitude  of  the  Rayleigh  wave  is 


just  given  by 


R = 2 


4.2  Implementation  of  Filter 


To  obtain  frequency  domain  information,  a DFT  with  a 256 
second  cosine  window  was  used.  As  explained  in  Appendix  B,  this 
window  prevents  leakage  of  the.'  signal  energy  among  the  harmonics 
of  the  transform.  The  window  was  moved  across  the  events  with 
shifts  of  15  seconds.  For  each  of  the  harmonics  from  64  to  19.  6 

seconds,  the  values  of  t-  , ft  , and  \ were  calculated  for 

every  position  of  the  window.  The  value  of  ellipticity  chosen  corre- 
sponded to  that  calculated  for  the  peak  amplitude  arrival  time  of  the 
Rayleigh  waves.  This  procedure  minimized  the  effect  of  interfering 
signals. 


4.  3 Results 


For  each  of  the  foxir  methods  of  calibration  outlined  in 


Section  3.2,  the  values  of  ellipticity  for  the  six  events  were 


calculated  as  outlined  in  the  last  section, 
cllipticity  magnitude  and  Figures  14  to  IV 


Figures  9 to  12  show 
show  cllipticity  angle 


for  each  of  the  calibration  procedures. 

Even  though  the  spread  in  azimuths  of  the  Aleutian  earth- 
quakes is  slightly  larger  than  that  for  the  Honshu  events,  the  ellip- 
ticity  magnitudes  arc  much  more  consistent  for  the  Aleutian  than 
for  the  Honshu  earth  qualms.  The  explanation  is  that  the  Aleutian 
earthquakes  with  azimuths  from  3.7°  to  8.6°  have  Rayleigh  waves 
which  come  in  almost  entirely  on  the  N/S  axis.  For  this  case,  if 
the  calibration  is  incorrect,  only  a scaling  error  occurs  in  the 
ellipticity  and  the  consistency  of  the  cllipticity  is  not  degraded. 
However,  the  Honshu  events  have  azimuths  from  35.5  to  39.7  . 
Thus,  each  axis  receives  a substantial  fraction  of  the  Rayleigh 
waves  energy.  In  this  case,  in  addition  to  scaling  errors,  errors 
occur  because  the  Rayleigh  and  I.ovc  waves  arc  not  separated 
properly.  An  error  in  the  ratio  of  the  x to  y seismometer  cali- 
bration looks  like  an  error  in  the  azimuth  of  earthquakes.  There- 
fore, the  best  estimation  of  ellipticity  magnitxide  will  come  from 
the  use  of  only  Aleutian  events.  However,  the  scatter  of  ellipticity 
values  for  the  six  events  in  Figure  10  gives  an  indication  of  the  per- 
formance of  ellipticity  filters. 
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Fig.  10:  Ellipticity  Magnitude  from  Average  Calibration 

by  Triple  Tuned  Model 
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Tlic  most  consistent  values  of  ellipticity  magnitude  (Figure  10) 
come  from  the  average  triple  tuned  model  calibration.  Table  II 
shows  the  estimation  of  ellipticity  based  cm  this  calibration  pro- 
cedure and  the  three  Aleutian  events.  For  periods  from  21  to  64 
seconds,  the  rms  error  is  less  than  6"‘<i.  Figure  13  shows  a com- 
parison of  our  results  with  those  of  Boore  and  Toksoz  (1969).  The 
difference  in  the  results  is  understandable  since  their  ellipticity  is 
for  the  Montana  l.ASA  station  and  ours  is  for  the  Kongsberg  VI. P 
station. 

The  ellipticity  magnitude  and  angle  is  more  consistent  for 
the  longer  periods  than  for  the  shorter  ones.  The  explanation  for 
this  is  twofold.  First,  the  effect  of  multipath  is  less  for  the  longer 
periods  than  for  the  shorter  ones.  At  the  long  periods,  there  is 
very  little  lateral  refraction  of  the  Ilayleigh  and  Love  waves  off  the 
ocean  continent  boundaries.  The  earthquake  azimuths  are  very  close 
to  their  great  circle  values.  Second,  the  calibration  signals  had 
more  energy  at  the  longer  periods  than  at  the  shorter  ones.  The 
increase  in  energy  improved  the  estimation  of  the  long  period  cali- 
brations. 

For  a perfectly  elastic  layered  medium,  the  theoretical  value 
of  ellipticity  is  90°.  Even  in  the  case  of  extreme  anelasticity,  the 
variation  from  90°  is  only  a few  degrees  as  determined  by  Boore 
and  Toksoz  (1969).  The  mean  of  the  ellipticity  angles  shown  in 
Figure  14  is  110°  and  that  in  Figure  15  is  115°.  From  those 
figures,  it  is  evident  that  although  the  fitting  of  the  tuned  circuit 


filter  gives  a reasonable  approximation  to  the  amplitude  response  of 
the  system,  it  introduces  phase  bias  errors  of  from  10°  to  15°.  If 

this  type  of  calibration  were  used  in  ellipticity  filters  for  mixed 
seismic  signals,  similar  types  of  errors  would  occur. 

llie  best  values  of  ellipticity  angles  comes  from  the  parabolic 
curve  fitting  calibration.  Figures  16  and  17  show  only  a few  degrees 
of  mean  error  for  all  but  the  shortest  periods.  Table  If  I shows  the 
mean  and  standard  derivation  of  the  ellipticity  angle  for  periods  from 
36.6  to  64  seconds.  In  general  for  the  four  longest  periods  shown 
m the  table,  the  fitting  to  the  average  response  is  slightly  better  than 
that  to  the  daily  response.  However,  for  the  very  short  periods,  the 
average  method  has  a systematic  biasing  error. 
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SECTION  III 


THE  ELEMENTARY  RA  YLEIC.H-LOVE  WAVE 
ELLIPTIC  IT  Y FILTER 


The  Elementary  Rayleigh-Love  Wave  Elliptieity  Filter 
is  one  of  the  simplest  Elliptieity  Filters.  The  filter 
operates  by  taking  the  outputs  of  a three-axis  seis- 
mometer station  and  generating  a Raylei gh-Love  wave 
decomposition  based  on  the  individual  frequency-time 
bin  data.  This  filter  has  been  described  previously, 
but  not  in  full  detail  which  is  given  here 
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1.  o 


Lnt  roduction 


The  Elementary  Rayleigh -Love  Wave  EllLpticity  Filter  is 
functionally  shown  in  Figure  1ft.  The  basis  of  the  decomposition  is 
that  the  observed  signal  is  comprised  of,  at  most,  a single  Rayleigh 
and  a single  Love  Wave  component.  No  allowance  is  made  for  the 
background  noise  nor  is  any  attempt  made  to  acquire  the  potential 
processing  gain  due  to  averaging  over  the  frequency -time  bins.  This 
filter  is  simply  a logical  extension  of  the  Polarization  and  Azimuthal 
Filters  with  the  important  addition  of  the  ellipticity  constraint. 


1.  The  inputs  from  the  seismometers  are  the  three-axis 

time  waveforms  )L[t)  , * and  corres- 

ponding to  the  E-W,  N-S,  and  vertical  instruments. 

2.  The  signals  are  sampled  and  digitized  and  the  descrete 
Fourier  transform  taken.  A time  aperture  of  about 
200  seconds  is  used  in  taking  the  DFT.  As  discussed 
elsewhere^,  this  duration  is  adequate  for  signal  periods 
of  as  long  as  60  seconds  if  a Cosine  I ader  is  used. 

The  frequency  components  of  the  Rayleigh  and  Love 
waves  present  are  expected  to  be  constant  over  this 
duraiion  and  to  have  almost  constant  azimuth  . 

3.  The  relationships  for  the  component  splitting  at  each 
frequency  are 

Xlf  = 

Y(r)  - 
zCf  - 


- £ If  R If  CoS  8 1 L C 

'R  l^) 


7 -l 


aiimlHaili 


' ■ ■ - -W*  •-'f  • - 1 


J-TT! 


in  which 

are  the  complex  Fourier  Transforms  of 
the  seismometer  signals. 

is  the  ellipticity,  defined  as  the  ratio  of 
the  horizontal  to  the  vertical  motion. 

is  the  complex  amplitude  of  the  Rayleigh 
event, 

is  the  complex  amplitude  of  the  Love 
event. 

B 9 „ 
' > /Z 

are  the  azimuths  of  the  Rayleigh  and  Love 
events. 

4.  The  solution  of  these  equations  gives  , 9,  , and  <5,  . 

Note  that  Rf^uj)  is  known  because  it  corresponds  to  the 
output  of  the  vertical  seismometer. 

5.  The  Love  wave  frequency  components  can  be  assembled 
together  and  the  inverse  DFT  taken.  This  will  give  the 
actual  time  waveform,  £Lt') . The  horizontal  Rayleigh 
Wave  component  can  be  similarly  determined. 


A similar  filter  can  be  constructed  in  the  time  domain  using 
the  analytic  function.  However,  there  appears  to  be  little  advantage 
in  doing  this  because 


1.  A double  DFT  is  required  to  give  the  analytic  function. 

2.  The  ellipticity  is  usually  specified  in  the  frequency  domain. 


The  instrumentation  of  the  Ellipticity  Filter  in  the  analytic  function 
(time)  domain  is  shown  in  Figure  19.  The  resulting  equations  would 


be  a set  of  equations 


the  Analytic  Function  or  Time  Domain 


^■0  = - £ it)  Kit)  Seng  - Lit)  co5  9Z 

V it)  = - £ it)  Kit)  ccS  <9,  +•  L[t)  Sen  £'^ 

^t)  = 7?£t) 

sirnilai  to  that  of  the  filter  in  the  frequency  domain,  but  £ it)  is 
harder  to  define.  lor  this  reason,  the  analytic  function  approach 
will  not  be  pursued  further  here. 


A mbiguities 
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2.  0 


As  expected,  this  elementary  decomposition  can  have  ambiguities, 
but  iu  actual  operation  these  ambiguities  can  be  removed,  for  the  most 
part,  by  using  additional  constraints  such  as 

1.  Knowledge  of  the  azimuth  of  one  component, 

2.  Behavior  of  the  reconstructed  wave  which  should 
have  the  required  dispersion  properties. 

The  potential  ambiguities  are  shown  in  Figure  20.  The  true 
azimuth  of  the  Love  Wave  is  always  obtained  as  well  as  the  true 
amplitude  of  the  Rayleigh  Wave.  However,  there  is  a spurious  Love 
wave  and  a spurious  Rayleigh  wave  which  can  be  present  in  the 
solution.  The  corresponding  amplitudes  are  tabulated  below. 


Condition 

Rayleigh  Wave 
Amplitude 

Love  Wave 
Amplitude 

Rayleigh 

Love 

T rue 

T rue 

K 

L 

Spurious 

T rue 

H 

E + ZER  sine 

T rue 

Spurious 

R 

— / 

/ — 

Spurious 

T rue 

R 

~ L ~ Z E R.  Sen  C 

As  noted  above,  these  spurious  waves  can  be  eliminated  when 
the  azimuth  of  the  Rayleigh  or  Love  wave  is  known.  Indeed,  given 


I 

I 

T 


S j is  the  azimuth  of  the  Rayleigh  wave 

is  the  azimuth  of  the  Love  wave 

X 

/f  is  the  Complex  Ellipticity  as  the  ratio  of  the 

horizontal  to  the  vertical  motion  for  the  Rayleigh 
wave. 


these  have  the  symmetrical  form 


+■  L Cod.  o4 


X - £ R Codo{) 

Y = ££  SiA/c<(  + L c.u^d% 


Then,  designating  real  and  imaginary  parts  by  the  subscripts  f and  jl. 
respectively 

~ +-/-r  Ccnzot^ 

X -u  = ^ Cod  o(  | + G>tJ 


+ /.r 

= (f  £^3>t  ( “h  j&Or'-* 

The  complex  rotation  can  than  be  used 


/l  = * r + 'j  ^ r - "h 


i ek^- 


V -Xi  + ^i 


= (£«):,  e^'+t.  £ 


(4) 

(5) 


in  which  and  V are  complex  numbers,  but  ^ 


£ . are  not. 


Then 
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Thus,  the  phase  of  is  U , while  the  magnitude 

ICE  H),  (£^J 


U k\  (W;J ' 


There  is  thus  a »■  ambiguity  in  oL^  for  which  the  sign  of  the 
amplitude  of  the  Love  wave  reverses. 


The  azimuth  of  the  Love  wave,  cL „ or  , cannot  be  solved 

A *- 


for  if 


i 

AM 

V 

i 

A 

i 

Li 

■ 

,(E*\ 

(£R,)‘ 

A* J 

This  situation  occurs  when 


The  Love  wave  amplitude  is  zero. 


In  this  case,  the  solution  is  meaningless 
anyway. 


The  Rayleigh  wave  amplitude  is  zero 


In  this  case,  the  vertical  seismometer 
output  is  zero  and  a simpler  set  of 
equations  can  be  used. 


The  Love  wave  component  (,'.)  is  in  phase  with 
the  horizontal  Rayleigh  wave  component  ( £ •s). 


In  this  case,  the  proportionality  constant 
( real) 


L = K (EfO 

has  a locus  of  solutions 

K '•  CoS  cos  °t- ^ X 

fd  ,o"v  cLf  + Sen  cC ^ Y 


These  are  triangle  relations  since  each  expression  involves 
only  three  vectors. 


A 


As  shown  above,  there  are  two  solutions  which  could  be  found  geo- 
metrically  by  erecting  JU.  & ° and  then  drawing  a circle  whose 


center  is  /)  and  radius 


(br\ 


One  of  the  simplist  relations  which  can  be  used  is 

A^r^(oL,-c<x)  - — J 

' ^ /k  m v ^ 


(£«)> 


If  o(  -oLx  — b is  a solution,  ff"  - Jj  is  also  a solution 


mm  -wAm^Mi 


«,  wwp 


There  are  thus  four  solutions 

2.  - ^ > c^l~^Z  = t>)  <^(_  c<^  -f-  b 

c>l  2 — T T 4 <3-  > C<  I - - Tf  — o(  I =~IT  - b + °tz 


when  written  out  are 


o4  , 

- cl  + b 

°6|  - 0-4  b + TT 

°^2 

- a- 

a +ir 

l = <2.  - b +■  TF  °i|"£L-b 

- CL-  = CZ+TT 

In  terms  of  the  original  angles 

<9,  = IK  -oi 

1 T 1 

®2  = ir  - ^ 

the  solutional  forms  are 

S,  = 3TT  - a-b 

2 

= if  - °~ 

e,  ^jL-a.  i-b 

z 

&z  - rr  - «- 


e,  ^Z-cL-b 
1 2 

ex=  - 
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Let 


c.z_b> 


b = ? - c 


d-TT-  a~; 


CL  =.  W - cL 


then 


<9,  = d + C 
Ox-cL 
6 1 ~ G%  ~ G 


0,  -d-hc. 

&3i  =.  cl  1-11* 

Q{-9„-  C -Tr  - ~rr  +C, 


<9,  ~ d-c 
Bx  - d 
dr6%-  -e 


9,  - d~c 

0%  =:  d -t-ir 
9,  ~ez  n -C  — 7T  ~ 7T  - C 


are  the  solutional  forms  for  the  azimuths. 


Note  that  if  ^ - *9(  = C — i s a solution.  (9^  - ^ — d C- 


is  a solution. 


It  remains  to  write  the  solutions  in  terms  of  the 


original  quantities. 
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1)  Relation  for  Q^_ 

Ox  - TT  - o( 


o«^  is  the  phase  of 


Xu  _ _\/ 

C ER)r  (ER)ju 


+ TT 

■b  O 


• . is  the  negative  of  the  phase  of 


m->  _ v; 

OLR\  (ER)i 


M'  _ V 

(£R)f  (£«)*,_ 

^*v  (e4 

i_Xr  _JU  ~j  +J 
(£R)r  f£R)ej 

k.- 


/>•  __  ^ 
(**)' 


■’■  ^ =.  —tala 


-/ 


Yu  -Yi 

C«)K  <ER)t 
XV  Jxl 
V«V  Ca?)v 


+-TT 

+0 


Quadrcnt  selector  is  done  by  the  signs  of  the  numerator  and  denominator 
Both  0j^  find  $ +TT  fire  solutions. 


5 -**H 


I 


(2)  Relation  for  <pj 


i! 


The  relation  will  bo  used 


,4-C/ru  (p<  | T.')  — 


-i4Z 


(EH), 


i! 


0 

si 


c<z  - 7T  - 

<*,  _ m _ (9,  -Ox) 


a^U'-oc^-  ^ ~ 69/  -<9i)j  ~ e— <fS 1 


Hence 


Gkl  69/  - ^ ^ ~ ^ 

(EK)r 


- — Cfh3  (P„  — Xf 


(e-X) 


r 


If  Q{~Q  -cL-  C.  i-s  a solution,  then  6t  — 0 ^ ~ C~  cL  i< 


solution. 


(3)  Relation  for  L_ 


I 

I 

I 

I 


T 


Mt. « 


4* 

«» 


it 


r t 


u 

11 


I! 


The  corresponding  equations  for  Ly-  and  L/u  can  be 
found  by  pairing  the  equations  (1),  (2)  and  (3),  (•!). 


Xr  Coock^  - CEfk)^  Gjw°<i  Cjk)g(z  4-  Lir  o<x 

lj,  y Ai^c<x  - tER)y  c< , o(2  4 ■ Lr  <u^xokx 


X f £<kv  4-  y s (£R\ 


or 


l~y  — Xh  G^olz  + ~ C£R)r  c^xlCoC,  -g(x) 


Since  ^ __  57 T __  @( 

1 Z 


<*/-<**  = JT  -CG, 


aiX  = T T - &Z 


A.  (r  ■=  - Xr  Cox?  Gx  4-  Cp  k 4^.  $;*,  — CER) ^ A~v\(Q  i -&*) 
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41  m. 


In  a similar  manner. 


X*.  db-doL^  - ^£f0o  (Vdo(j  G<X)o(^  -f-  LL  Ccv 

— L|  a.  jdliw  o(^  ~ (^  E Ej)j^  ’&L*r\  °(  j jd-t-m.  •/-  In  a>  £■*-'>**?' at 


C<A 
.2. 


X,o  Q.-tx>o(  ■+  .u  i,  di*  o(  - C^L  (c^l  - )+J-~X> 


(7) 


i-I-  " Xxl  4^o<*  CoxLiot,  -o{:,) 

or 

,L.  A->  — — X-O  C.(hi3  ® 2.  *^“  ^ ^ ^ — (j£  f?)  ^ AlA-r^  (Q  { ~Q 


The  final  relations  for  ^ and  /.t  show  that  if  /.  is  the 
solution  for  6^  = C , Q ^ - cL  , then  — L is  not  the  solution  for 
o,  = C + 7T  , @2  - ci  + If  because  although  i9,  and  6^  change 
by  7 T , the  difference  does  not.  The  two  relations  for  /.^  and 
given  above  can  be  summarized  as  follows: 

L * ~ \ cos  Bz  + V 5 in  0X  ~ [ £R ) sin  (9,  - Q^) 


<)l 


s-T”  , prVvj 


(4)  The  pattern  of  the  solution  can  be  described  as  follows: 


Fo  r 


0 1 = d--t  C (true  Rayleigh)  _ )<  Co-d.  cL~\t  / _ £/?  C 

(true  Love) 


This  is  taken  as  the  true  relation 


Si  -Qx  = e. 


and  will  be  designated  as  L £. 


L - 


For 

&/  — ci-C  (Spurious  Rayleigh) 

Qx  - d~~  (True  Love) 

Q,-0X  = -C 

L—  x c*-o  t£j- V zL -t  £<? 
- i-t  + ^ c. 


For 

(9/  tfUC  (True  Rayleigh) 
$2  _ +TT  (Spurious  Love) 

#1  - <9 Xrz.  C 


£_  =.  •+•  x Ccm  d.  — y dL-t- 

- - Lt 

For 

<9i  = <^-e  (Spurious  Rayleigh) 

(9x  = cL-+Tr  (Spurious  Love) 

(9/  -^  = -G-tt 

i_  - x C-Ots  <sC  — y d-  — E.Rj  At^*y.  CL. 

U--L.il  — 2.E.R.  A^-Orx.  C. 


92 


Numerical  Example 


I 


Select  ££  = * + ^ 

L , \+*i 

9 ( b20° 

= 30° 

Then  X = — H R.  £lL^,  Q \ _ i_  C^si-. 

= ^ — /QjUw  20°  - cos  30°  ^ 4-  20  — X C^l,  30  ^ 

- - 1.  2080  - ^ 2.  0741 

y = — E O^crQ  0 1 + L <£^k,  £9  X 

- C~  C'Crd-'  20  + sin  30  ^ C ' cos  ^ 30°^ 

= - 0.  4397  +-i  0.  0603 


ir-  _ 

(EOr 

— 

(e^V 


_y±_ 

(er): 

Xi 

o^x 


= - tan-'  - °-  4397  - °-  0603 


- 1.  2080  + 2.  0741 


- ian~  (-  0.  5T73) 


Fourth  Quadrant 


Qn  = 30°  and  210° 
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os  C $(-$*):= ~~  X-r  _ 

<^JT 

or  fi*  = 30° 

cos  (Si  -Ox)  - 0.  4397  cos  30°  -!  1.  2080  sin  30c 
= 0.  9848 

= 10.  oo°  , S, -Qx±  io° 


10.  00 


ror  #2.  = 210  , (180  + 30°) 

(9,  _ Q*.  = 170°  ^ - 170° 

and  - 170° 

9\  = 20°  0,  = 40° 

= 30°  02.  - 30° 

$1-0*=  - 10°  Qi-Qx  s + 10° 

(9,  = 40°  01  = 20° 

9*  - 210°  =:  210° 

= - 1-70°  $|-5^=  “ 190° 

The  Rayleigh  component  is  at  20°  (true  angle)  or  40°( spurious  angle). 
The  Love  component  is  at  30°(- rue  angle)  or  210°  (180°off). 


rjuJcLi 


i v;,. 


Ah  Ir  cos  Ox  + tyy  Kin  0 z “f£ft)^sin  (Q\~9*) 

A ^ "—Xx-cos  9%  **  sin  Cdd')  ^ sin  C@t  - & z) 

9t  ~ 20°  A \r  - \.  2080  cos  30°-  0.  4397  sin  30°-  sin  (-10°) 

&%  = 30°  =1,  0000 

#,  -0%  - -io° 

A = 2.  0741  c.os  30  -I  0.  0603  sin  30°-  sin  (-10°) 


Q(  = 40° 

n Air  s 1.  000  - 2 sin  10°  - 0.  6527 

^ 30° 

Air  = 2.  000  - 2 sin  10°  = 1.  6527 

+ 10° 

40;)  A = _ o.  6527 

= 210°=  180°+  30°  Air  = _ i.  6527 
- - 170°=  180°+  10° 

$/  = 20°  Ah  “ - 1.  0000 

@x  =•  21°° - 180°-l  30°  Ah  = - 2.  000 
8l-®X=  - 1 90°-  + 180°-  10° 
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SUitaJ^ 


Spurious  True 


40! 


30" 


0.  6527  + 3.  1.  6527 


True  Spurious 

Spurious  , Spurious 


20° 

40° 


I 

i 


'210 


i -0.  6527  - ^ !•  6527 


Asa  check 

$|  = 20° 

X 

= ( - sin  20°  - cos  30°)  + j(-  sin  20°  - 2 cos  30°) 

~ 30° 

= - 1.  2080  2.  0741 

y 

- ( “ cos  20°-5  si«  30°)  _ cos  20°4  2 sin  30°) 

L = 

- 0.  4397  +j  0.  0603 

0,  - 40° 

X 

- (-  sin  40(/-  0.  6527  cos  30°)  4'j(-sin  40°  - ],6527cos  30°) 

(9Z  = 30° 

= - 1.  2080  2.  0741 

£R  = J+i 

y 

~ (-cos  40  4 . 6527  sin  30)  l^(-cos  40°  I-  1,  6527  sin  30) 

L 0.  6527  +^' 

1.6527  = -0.4397  +^  0.0603 

(9,  - 20° 

X 

- (-sin  20°  + cos  210°  ) 4^  (-sin  20°  1 2 cos  210°) 

<9^  = 210° 

= - 1.  2080  -i  2.  0741 

£&=!+% 

y 

= (-  cos  20  -jsin  210)  4-^  (-  cos  20-2  sin  210) 

L 

= - 0.  4397  4-j.  0.  0603 

96 


.wnaini  


SECTION  IV 


THE  ELEMENTARY  RA  YLEIGH-RA YLEIGH  WAVE 
ELLIPTIC  IT  Y FILTER 


The  Elementary  Rayleigh- Rayleigh  Wave  Ellipticity  Filter 
is  a simple  Ellipticity  Filter  which  is  based  on  using  the 
three-axis  outputs  of  a seismometer  station  to  generate 
a Rayleigh-Rayleigh  wave  decomposition.  This  filter  has 
been  previously  described,  but  not  in  the  full  detail  given 
here. 


1.  0 Introduction 

The  Elementary  Rayleigh-Rayleigli  Wave  Ellipticity  Filter  is 
functionally  similar  to  the  Rayleigh-Love  Wave  Filter  except  modified 
decomposition  equations  are  used.  The  formulation  is  best  done  in 
the  frequency  domain  as  described  here.  No  allowance  is  made  for 
background  noise. 
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2 . 0 Performance  When  a Single  Rayleigh  Wave  is  Present 


When  a single  Rayleigh  Wave  is  present,  the  received  signals 


X n - ft  Stn  0 


V,  = X CoS  0 


ft 


The  azimuth,  9 , and  the  complex  amplitude,  & , are  uniquely  obtain- 
able from  the  observed  signals 


ft  = 2 

xo£v<  0 3 J Z 


cos 


b = y0  / z 


'l 


1.  1 


! 1 


1 

f ) 

I 


* 

i 


The  presence  of  a single  Rayleigh  wave  can  be  distinguished 
from  two  Rayleigh  waves,  although  it  may  not  always  be  distinguished 
from  more  than  two  Rayleigh  waves.  The  test  can  be  made  by  exam- 
ining 


r - x.*  f V/  - ** 

and,  if  T = 0,  deciding  that  one  Rayleigh  wave  is  present. 

The  basis  of  this  test  is  that,  for  two  Rayleigh  wave  components. 





x,  = 2, 

so»  Q ( t 

S ih  6Z 

yo  ■ R, 

CoS  +- 

^a 

cos  Gz 

z = ft' 

T = X * 

o 

r V - 

w 2 

X 

= 2 K,  R z [ SiK  @t  si,h  &2  + cos  coS  &2  ~ tj 
= tek  - 9X  ) - ij 

Hence,  if  ft,  ^ 0 or  / 0,  the  test  quantity,  T*  , is  zero  i 

only  if,  &{  ~ &x. 
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, and 


3 . 0 Performance  When  Two  Rayleigh  Waves  are  Present 

The  direct  method  of  decomposing  two  Rayleigh  waves  has 
been  previously  presented*.  The  main  results  are  reviewed  here 
with  additional  comments. 

When  two  Rayleigh  waves  are  present,  the  basic  relations  are 

X = sin  0.  +-  H . Sen  B~ 

O ll  < 'L 

CoS  0(  4-  cos  &2 

* « *,  * 


in  which  the  ellipticity  factor,  assumed  known,  is  included  within 
the  observables  Xe  and  Ye  . In  terms  of  the  real  and  imaginary 
parts,  these  equations  become 


X0r  - 

K 

Sen 

0, 

+ 

Hzy  Sin  B2 

^ ■ 

K,i 

Sch 

t 

K ■ sin  0 . 

* 

Multiplying  the  first  equation  by  ^ • and  the  second  by  %iy  and 
subtracting  gives 


= ( %,L  - *ZL  sin  0. 


Defining 


P, 
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gives  the  compact  form 


(*.;  - 

P,Xo>' 

) = 

^ ^ ^ SLn  ^2 

In  a similar  manner. 

the 

y. 

equations  can  be  paired 

t. 

= 

K5„ 

Cos 

^ cos  a* 

V. 

= 

Ru 

cos 

B t t K^i  cos  v9^ 

to  give 

^ - 

V 

1 oy 

) . 

The  relations 

- 

*, 

- 

% 

1 ►- 

-S' 

/S  . 

L 

= 

- 

can  be  used  to  reduce 


% ti 

fc... 


(*i  - *0  ' - b,k> 

**  <r 


'R’ll 

lv 


Hence, 


c%°-  - p,  *„P 

( V„;  - P,  V„  „ ) 


i.  %■  l H , y ) SUM  Qz 

( ~ Pj  £y  ) CoS  9, 


Squaring  and  adding  then  gives  the  quadratic  form  for  ~Pf 
dP(^  +■  b P i +-  C = £> 

By  a similar  derive. tion,  the  same  quadratic  equation  is  obtained  for  p, 

'P  = p • / p 

-2  r2c  ' r2r 


Thus,  P, 

and 

Pj  are  the 

roots  of  the  quadratic  equation  in 

a 

= 

2 

+ V z 

f or 

- ^ 2 
o>- 

b 

■= 

* L 

Ki  “ 

^ V„;  - V J 

c 

= 

xo* 

r V * 

oo 

2, 

Aj  . 

o» 

Note  that 

a 

- 

<X 

Sun  9 

+ V 

5Jh  0)  ^ Cr,it  . cc; 

- 

■2*V 

E.,  i 

- «•)  + 0 

+ p ) 
ir  l^z^J 


t = 2 [(%  l(,  t - CP(),  cos  9 + cos  cosfif  + K^cosS, 


(X)k,  Sun  9,  +•  Sen  9*  X &(l;  Suh  <5,  t t*  Sih  6^) 


~?r-  "*i  ^?.  \.  f;  * • 'V- v^T"'  'T'  UL1/. Wi.*,  *)K‘ 


= 2 LX 

f 1 

- .2 
— ros 

— s«.’h*  fi,') 

r 1 

- c 05  6( 

C^S 

“*  5ch  *5^  Sc<) 

®2  ) 

Ci 

- cos  <9^ 

COS 

:^j  S c h S 6 

n (9( 

*«■ 

*,i 

■ (l 

- So  n 2 ) J 

•>? 

1 

II 

*ai 

+ K,.  R 

,L  )f  Cl3S  ^5i  ~ 

*0 

= - zx>,i  k^l  C c<>s  (e,  - sz)  - 0 

Hence,  if  0,  = 6Z  , a.  = b -■•  C is  zero.  This  degenerate  case 
can  be  avoided  by  first  testing  for 


T 


X. 


and,  if  T = 0,  using  the  single  Rayleigh  decomposition. 
With  P and  P solved  for,  6,  and  can  be  found  from 

X„; 


totv,  6 . 


? X 

' I 


<K 


- r, 


and 


ta.W  6 


- p,  *= 


>U  - ?*it 


oe 


These  have  a 7f  ambiguity. 


If  9t  = 9{  , then  0,  = 9f  f TT  is  also  a solution,  and  if  ^ 

then  9^  = 9^  t 'if  is  apparently  a solution. 

These  ambiguities  in  azimuth  can  be  resolved  by  using  two  of 
the  three  relations 


cos  9a 


(-  ^ ) 

- ?,  Re  L-  *ll) 

3,  ( z~) 

r,  Re  Cz) 

- P,^£L~V/e) 

3*12) 

- ?,  2e  C 

3,  t- %) 

- P,  Rc  (ry/e  ) 

3,9%) 

“ ?.  6" 

or  testing  the  quantities 


m^F**!*-  ^^mrr  m^-^.^rSTT^I.V^b  ■i^4*i'? -^AW* .-V-^i J UIP^.^IW 1 1 .1 , |. -^«J  -",,1'A.l^U  V • I ! J P^ILIIJHI | 


vam*:  fjms93pm*T iw» 


and  assigning  the  quadrant  by  the  rule 


B 


+ 

+ 

+ 

+ 


+ 

+ 


+ 

+ 


Sine 

+ 

+ 


+ 

+ 


Cosine 


+ 

+ 


Quadrant 

I 

II 

m 

IV 

IV 

m 

n 

1 


The  amplitude  and  phase  of  the  Rayleigh  components  can  be 
resolved  as  follows. 


If  now 


% 


/file 


iRt'l  Cos  + jlzjsuhfy  = ^,K  + -j^ic 

Z z - l&J  * = IR,J  cos  ^ i-  j I sui  tfz  - Rjy  r 


,4 


it  follows  that 


ta.fr  (j) ] 

tam 


F z 
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except  that  there  is  a IT  ambiguity  so  that  if. 


4'  0 , 


d>  = + TT  is  a solution 

Y / T / 


<pi  " & > ^ 


is  a solution. 


To  resolve  these  ambiguities,  the  magnitudes  of  A?,  and  H are 


solved  for. 


The  magnitudes  Z and  Z can  be  found  from 

Iw  2m 


Zt  + Zz 


z,  ei<f'  + z 


£ -h  z e 


it*,  -^) 


+ ^ \rn^°S  ($!-<?,) 


The  imaginary  parts  of  this  equation  give 


Zlcos  fi  - &y  si»<ft  - Zzm  SL«($Z  ~(p,) 


Zfc  CoS  P ! - y,  si*  4‘l 


WWWr» 


If  K is  to  be  positive,  $ must  be  chosen  as  fz  or  ^ +■  7 T 
do  so.  This  selection  is  unique. 

Simila  rly 


% 


I r* 


X 


■t 


XT 


•Si^o 


f,  - 


CCS 


to 


requires  the  same  selection  between  <p(  and  tfJt  * 'll 

What  may  not  be  apparent  from  this  direct  method  is  that  the 
solution  is  ambiguous  if 

^ It  / ^Jiv  = / zv 

or 

f>  = P 

Then 

= *,„e^  = n^ei'e 

**  * - -R,„ei<P 

since 

f,  = f 

Then, 

XD  = 0,  + ^ 

<U  * Co*  + ^a»n  CoS  ^ e 


- — IIHMII  mil—  WMlWii^iiSfriW^iWiigii^y  .i^jpwk^^, 


or  all  the  observables  have  the  same  phase  angle.  While  Q can  be 
solved  for,  the  equations  are  degenerate 


Xo*  = 

Siin 

+ 

&z 

y0„  - 

cos  0, 

+■ 

cos 

8z 

a£  = 

OK% 

+■  1 

> 

There  are  three  equations  for  four  unknowns. 
Z si  2 w Z 


if  x„  * + y/6  - 2 

Cry\  0*v\  D W) 


is  formed 


X Z -h  y z 

or*  r 'a*. 


- % 


= Z *2  2, 

!►*>  /J»l 


£ 


COS  ^ <9,  - Gj,  ) ~ ij  ^ O 


and  one  can  select  any  0,  - 9^  as  a solution. 

This  case  is  quite  unique  and,  thus,  not  too  likely  to  occur.  To 
assure  ourselves  that  there  are  no  other  ambiguous  solutions, 
indirect  method,  discussed  in  the  next  section,  is  used. 


the 


vii’.*  Tw  ■ ‘ ■ 7 - ' Ys  iTTr  w ™ ~r i-'wt' y»7 v*-^'.-  • . 


4.0  Performance  When  Two  Rayleigh  Waves  are  Present: 
Indirect  Method 


Additional  insight  into  the  uniqueness  of  the  solution  can  be 
obtained  by  a somewhat  different  approach.  The  direct  solution  has 
a jump  discontinuity  as  9 ->  <9,  . To  avoid  this,  let 


v = - fcj 


lyi*.  i-  V )/  Z 

i o.  -'•')/  z 


since  these  variables  have  no  such  discontinuities. 


Then 


^/z  L Si-h  P(  + 5ih  +■  2 L Cos'  — Sit!  6^] 

'UU  [ CoS  fi,  4 COS  9Z]  -f  V/zLC"S  >9,  ~ COS  <9j. 


The  basic  relations 


Sun  + SC*\ 


, a,  4-©^  [b,  ~ez\ 

Z s.h(  7 J cos  ( J 


~~  5u*"i  & ^ ~~  £ Suh 


««  (M--) 


2 »s  ~ (^-p-) 


CoS  6 ( - Cos  0; 


t ^3 


suggest  another  change  in  variable 





oi 


9.  -r  9. 

1 A 


9,  - 9X 


r * 

This  gives  the  set  of  relations 


— 

z 

Son  rd  COS  yS 

V'  C05d(  Sin  jl 

Y„ 

= 

2 

cos*:  cos  /3  — 

l/  5th  5th  3 
/ 

In  terms  of 

the 

real 

and  imagina  ry 

parts 

- 

5tn  «:  cor  p>  + 

CoS  5th  p, 

3 

coS  K.  cos  — 

J y Sin  a.  Sin  3 

*-  + i 

- Xy  COS  f!>  ( 5 in  ot.  i 

JC5Sc<)  + V y.  Sin/i  Ceos  a.  - ^ sih  «*'.) 

e 

'iu 

cos  cC  - 

^ So  h ch 

i e 

-i* 

= Sen  pt  +■ 

~4  CoS  oi 

O 1 i0 


'Zy  CoS  p ^<9.  -f  Vy  Sin  p>  e 


- v|  oC 


In  a similar  fashion 


X„  v 

Oi. 


i ioi 


- ce>s 


e 


- (oi* 


- -/<=£ 


+ V;  sih 

Aa 


/* 


e 
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- — -- 


rhen, 


^oi  + ^ ov 


/3  6? 


^OV-  + i 't  L 


Xoi  + 1 ^OL 


cL  or  oL  + V~  is 


the  phase  of 


>w  + riV, 


+ 4 ^oi  A 


‘OL  r A 


Hence, 


ta.^  cL 


Once  oi  is  known,  the  other  quantities  v and  j$  can  be  solved  for. 


The  solution  breaks  down 


If  y6  =0.  This  corresponds  to  or  the  cast 

in  which  only  one  Rayleigh  wave  is  present 


In  this  case 


Vi/  , 

' *1 
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This  corresponds  to 


Zir  + %3.r 


% ti  ~ tZ&L 


%,L  + '&zl 


or 


"l^iL  ! 
/ 


Physically,  this  corresponds  to  the  two  Rayleigh  waves 
having  the  same  phasing. 


These  are  the  only  possible  ambiguous  cases  that  can  be  present. 
Numerical  Example 


■ ' + zi  = 

2. 23&I  @ 

= a**;  - 

5.  oooo  (o' 

S3.  l3o  1° 

*1 

0 

o 

II 

II 

a 

o 

£ 

■ i 

= (_  I +■  J2  j ) 5th  20°  + (Z  +•  4-j  ) Sth  3o° 

= C Sim  H 0°  +■  3 Sin  3o°)  -f  £ C 2 sin  20  + 'f’ 5th  3o°) 
- I.  $420  t ^ 2.  (,$40 


C I 4 2j)  cos  2o°  + iz  +•  4j)  cos  3o° 

C cos  2o°  + 3 cos  3o°)  +-  ^ ( 2 CoS  2o°  4 4~  cos  Jo°) 

3.53  7?  + 4 S'.  3+3 T 


y 


(3.5378 r 1 ( 1.  8420p 


(4) 


a, 

= - 0. 09101 

b = 2[  (4)  (6)  - (3-  5378)  (5.  3435)  - ( 1 . 8400)  (2.  6840)  ] 

- 0.30368 

C * (5.3435)2  4 (2.6840)2  - (6)2 

= - 0.24315 


aP2  + bP  tc  =o 


has  the  solutions 


"P,  = 2.002, 

(2.000) 


Pa  = 1.3339  (checks) 

(1.3333) 


s 


- Tjge(-T') 

L-tt)  - p;  R=(-f). 


tah, 


-I 


- Z.Ll^O  +•  £ 2 .Ooz)0-  $+Zo) 
-57  3435’  4-  ^2-ooa)(3- iT3  7^) 


— taw. 


_/ 


f-l.ooo  0 

+ I ■ 732/ 


5 taw,"'  0.57733 
3 ° and  J2/o° 


However, 


= 1.0000 
3 = 


1.73  21 


c 


since 


iU-V)  - P,  Re(-*-) 


h 

B 

C 


- 5T.  3*35 
/•  7311 

is  + 
is  + 
is  + 


- Cz.oo2.X~  3,^37?) 


only  30.  00  is  the  correct  solution 


However, 


G. 


~ t, 


-i 


ah< 


Since 


(~2..t>8+o)  ~ ( /-3333)^-/. 


~ 3.  3435  - 0.  3323)(~  3.S3*fO 


— taM, 


-O.  22#ol 


H 

B 


- o.c>zt> 

to.h.”^  31>Z^7 

2o.oo0  and  2.0a 

- o.  22  go  I 

- O.  6264-5' 


c = u(~r) 


- sr.  34  3 s - 

~ o.  6 2 6 * 


- P,  Re(-i) 

Cl . 333&)(-  s-terO 


A is 
0 is 
C is 


only  20.  00  is  the  correct  solution 


1 16 


Hj 

5 .^1 


4 


amplitude  ftf  = 


iLM  - ?x  Ke(z) 


for  * (,3.4.3+° 


SUv  frfi  - CoS  ff, 

6 - 0‘ 3333^6^) 


‘ j 


amplitude  ft  - 


Sin,b3.43+°  - 0.3333>  Cos  (,3-43+c 
2-2347  (checks) 

4 - C\.  3333X4) 


sin,  243.+3+0  - { 1 . 33 33)  cos  243.-4-34' 


for  ft  - 243.4 34° 


1 

a 


= - 2-2343  (checks) 


Sign  of  /)j  reverses  so  that  this  is  one  solution 


i 


ft,  = 2.  2347  ZJL£_!£±i  (checks) 


: 


amplitude  /)2  = 


4 - 2.poooC  +) 


Sin  S3.  I3o°  — C2..0000)  Cos  5"3./3ijc 


for  (f>  = S3  ./3o° 


S.0000 

S. 00 no  / ^3-  <3o°  (checks) 


1 


Hence,  the  unique  and  correct  solution  is 


e, 

0, 


= 20 


'X 

ft, 

*2 


= 30 


o 


= 2.2361 


4349 


= 5.0000 


Lh 
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SECTION  V 


AN  ELLIPTICITY  FILTER  CONFIGURATION 
BASED  ON 

MINIMUM  NOISE  MAGNITUDE  ESTIMATORS 


It  is  sometimes  desired  to  estimate  the  magnitude  of 
an  event  from  surface  wave  data,  when  the  surface  waves 
are  partially  masked  either  by  the  main  surface  wave 
arrivals  or  by  the  late-arriving  Rayleigh  wave  coda  of 
another  seismic  event,  and  when  noise  may  also  be  pre- 
sent in  the  record.  A minimum  noise  solution  for 
magnitude  is  developed,  based  on  an  Ellipticity  Filter 
processor  in  which  a flexible  three-component  signal 
waveform  is  derived  from  surface  wave  data  from  pre- 
vious signals  from  nearby  epicenters.  In  the  case  of 
interfering  Rayleigh  and  Love  waves  from  a known  epi- 
center, the  solution  for  signal  magnitude  has  been 
derived  in  closed  form.  In  the  case  of  interfering  Ray- 
leigh wave  coda  from  unknown  directions,  the  solution 
is  obtained  by  a single- variable  search.  It  is  expected 
that  these  solutions  are  considerably  more  robust  than 
those  of  the  Elementary  Ellipticity  Filters  predicated  on 
low-noise  data. 


1.  0 Introduction 

The  problem  is  to  estimate  the  magnitude  of  an  event,  based 
on  the  Rayleigh  and  Love  wave  signal,  when  the  event  is  mixed  with 
the  main  surface  wave  or  the  late -arriving  surface  wave  coda  from 
another  siecmic  event.  The  approach  is  to  find  that  Ellipticity  Filter 
decomposition  of  the  signal  and  the  interference  which  minimizes  the 
residual  noise  in  the  records. 

The  method  utilizes  the  assumption  that  the  actual  signal  will 
differ  deterministically  from  a "prototype"  signal  (empirically  modelled 
from  similar  events  from  the  same  region)  as  a result  of  differing 
origin  time,  range,  azimuth,  and  Rayleigh/Love  energy  ratio  as  well 
as  magnitude.  The  method  further  assumes  that  the  interference  is 
either  a combination  of  Love  and  Rayleigh  waves  arriving  from  the 
epicenter  of  the  other  seismic  event,  or  that  it  is  late-arriving  coda 
consisting  only  of  Rayleigh  waves  from  unknown  azimuths  corresponding 


to  various  multipaths. 


.w****8W9WS«n^^ 


2.  0 Derivation 


The  three  observables  are  the  outputs  of  the  three  x-axis 
seismometers.  As  before,  the  solution  is  obtained  in  the  frequency 
domain.  Thus, 


^ j^f j % are  the  Fourier  transforms  of  the  time-domain 

observables  corr esponding  to  eastward,  northward, 
and  upward  motion  being  positive. 


The  decomposition  is  based  on  the  prototype  waves  being 

for  the  amplitude  of  the  prototype  Rayleigh  Wave 

L. p for  the  amplitude  of  the  prototype  Love  Wave 

for  the  amplitude  of  the  interfering  Rayleigh  Wave 

/-£  for  the  amplitude  of  the  interfering  Love  Wave. 

These  quantities  "R^  , , R t , L£  , X , j , and  ^ are 

all  complex.  The  notation  will  be  used  that,  if  Q,  is  complex,  Q.*" 
is  the  complex  conjugate  and 

ft  ' %r  + ‘?i 

For  the  actual  event,  the  azimuth  associated  with  each  wave  is  desig- 
nated by  0 , 9^p  , , and  8j;t  respectively.  The  magnification 

of  the  prototype  Rayleigh  and  Love  waves  are  designated  by  f M 
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7'  -"j’-.rf  ■'  ■ . - 7-^7 crr  — t-tt t tv’  */  n '■»  j i-Ti? »■■ 


The  noise  is  denoted  by  N x-  , N,  , and  N corresponding  to  the 

<i  o 

three-axes  of  the  seismometer. 

The  range  dependancy  effects  are  taker-  care  of  by  the  inclusion 
of  an  exp( A.  ) factor  in  which 


t'i  t */V) 


being  the  frequency  of  the  Fourier  transform  bin, 
being  the  delay  in  origin  time, 

being  the  increase  in  range  from  the  prototype  reference, 
being  the  phase  velocity. 


If  no  subscript  appears  on  A , it  is  implied  that  A and  V corres- 
pond to  the  wave  type  associated  with  the  exp(-^A)  factor. 

The  decomposition  equations  are  then 


X s ftp  & Sin  Gyp  f-  'Rj^.  SCm  6 y^) 


yp  p yp  1 x yx  J 

- (Mjp  e r COS 9 + tx  c*6tt)  + 


t?,p  e 


CC5  Qrp  + nx  cos  Qyr) 


+•  Lp  € ^ Sin  Qjp  +■  l~x  S-'n  d£x)  + 


m..  -r,  e 

rp  P 


+ Ei  + 


••VH/n/Si&itftffrpf  iaKSr^ttJi 


The  decomposition  is  based  on  minimizing  tlie  noise  sum  expression 

N2  = f NkN*  -h  * 

by  selecting  Ny  , , %£  , l~£  , 9r£  , and  S^£  to  do  so. 

The  factor  VJZ  is  permitted  to  allow  for  the  situation  in  which  the 
noise  output  of  the  vertical  seismometers  is  substantially  different 
from  that  of  the  horizontals. 

The  solution  is  formally  obtained  by  writing  the  noise  sum 
expression  and  setting  the  derivatives  equal  to  zero.  To  perform 
the  algebra  with  minimum  handling  of  subscripts,  a change  in  notation 
will  be  used. 


ft 

' Kr 

e»  • 

c 

= *> 

<5=  - 

Mc  ■ ■'V 

J> 

* 

&S>f> 

The  decomposition  equat  ons  are  then 

X = - £ ( £ si-h  Of,  + f\  si> 7 & # 3 

“ "6  cos  C cos  Bc  (2.  + N x (2a) 

V = - E [ M?  j)  £ * sin  8- j,  +-  8 cos  8^ 

_ i A 

t B sin  Sg  + Mc  C.  Sun  Gc  6?  + (2b) 

Z « j>  e~LA  + a + *iz  (2c) 
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. . /P 

Ny  Niy 


X X t £ J>  J>  sihZ  8 M * + £ E*  f\  f\  *sm* 


J>  D 

+ 8 B*  Cos*  eB  + C C^ces*  ec  MXC 


6. 


* — e. , 


cos 


+ Z He  ^ X E * A*  si*  + X 8 
+ Eh  si^Q  co£  <9  + c cos  6r  M e~LAL  X * 

n H c c 


+■  e+  S^h  0JJ  t E*/}*Sfch  + B*  CoS  &BJ 


+ E 3>  e*t4  si"^  [x*  t A*E*st*&t  + B*cos8bJJ 

(3a) 

V 

same  as  Nx  N*  , except  that  X becomes  V , 

sin  becomes  cos,  and  cos  becomes  -sin.  (3b) 


N»< 


* ^ ^ A*  h 2”  Z * -h  2 Re  £-  ZJ)*g  M 


13> 


y -LA  %■  -i 

+ J>  e /I  M J (3c) 


2 

The  derivative  of  N with  respect  to  is  as  follows: 


a/v; 

bM. 


W*k3>P*tfp  + 2«e{j>  e'tA(7)*  - z*)j J 
+ [z  £ E*.P  sim^dp  Mp  + 2£eJj>  e E si h&p  (X 

i~  ft*  E.*  sin  6 + B^CcS&a  + C*e+c/i  COS  &c  WC)]J 


* 


f [*£  £*;>  V*MX&p  MP  + ^[^e“  EcosG^d 


*' 


•f-  fi*E*cos8A  - B*sL*9b  - c*e^^sCh&c  Mc)]] 


1 Z 3 


MAMMiiiiiliMii 


atiiliii-'li  ITI I 


This  can  be  reduced  by  noting 

- JP^Cw>tE£V 

d Mp 

+ ZHe^E  J>  e'L*  C*  efL^  (sin  cos  &c  - cos  O^sinQ^ j 
t 2 Ze(}  e’LA[EX*sii  + Ei^cosB^ 

+ A*E*E  (sin  6 sin  8^  + cos  Op  cosB^) 

t £ B*  ( S£n  fijj  CoS  9g  — COS  sin  &■  g ) 

t ua(/)*-  2')]? 


^ M, 


2 P J>*  (U2  +■  E E*)b4p  +■  2%e(EPe  LCi C^e+tA) sih^-^ 

t-  2fc?[;P  £X*sin^  + El*  cos  Bp  t A*  E*  E cos  ($^-6 

+•  £B *sin(Q  -e^  -h  W*(A*-Z*)jj 


j> 

(4) 


Note  that  the  first  term  on  the  right-hand  side  of  Eq.  (4)  is 
linear  in  and  has  a known  coefficient;  the  second  is  linear  in  A) c 

and  also  has  a known  coefficient  (which  is  very  small  since  <9p  fc  ); 
the  third  is  dependent  on  the  unknowns  related  to  the  interference. 
Now, 


brt, 


2C  C*tAc  cosz  ec  + 2/?e£c  e~L*Cos9c  [x* 

* E*  J)*  e+l£t  sin  6p  Hp  -h  E*fl*sin  &A  + B*CosOr 


]} 


+ 2 C c*  Mc  sin  X 9C  + 2 £e  C e Sch  Bc  [ Y * 

4-  eW^os  &.M.  + eVc*®*  - bVm 9e:( 
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ii 


n 


o 

0 


b H 


ra 


^M, 


«c'flc  t 2 Re  [Ce  J>  e“ae,e’f‘astr,/'ej)  - 0C)M 


+ C e-L*  [ X*  ces  ec  - Y**a8c  t E*fl*sU  (B„-8C') 

- e<  (5) 


If 

t B cos 


Equation  (5)  has  the  same  properties  of  linearity  as  Eq.  (4), 
with  the  roles  of  Mc  and  reversed.  If  Mc  = , 

~bNZ 


r al  ^ 


~bn 


M 


M, 


2>,V: 


bfl 


= 2[j>  ;p*(U/*  t £Ef)+CC 


+ 2Zz  (£*£  C e"tA)Sth^  Bp  - <9C)J  M 


f 2 Re  ^Remaining  terms  of  Eqns.  4 and  5 which  do  not 
^ include  M~) 


{6) 


In  a similar  manner,  the  partial  derivatives  of  N with  respect 


to  A , Ti  , 0^  , and  may  be  found.  Since  the  partial  derivatives 


of  NZ  with  respect  to  the  real  and  imaginary  parts  of  H ( and  <xt 
respectively)  are  both  real,  it  is  equally  useful,  but  more  convenient, 


1 

j 


to  find 


T £ 


^ o*v 


'la,  l 


s zu*L/\ 

+■  2 [ E E * /)  sih  * 9 

f\ 


Z + JD  D 


^*0*  t X E*scn  8h  + 5 i/i  S c0s/96  £* 

Sen  Mjj  £ *sih  6 ^ 


+•  C cos  9C  Mc  e"tA  £ %ih  eh  + 


+ z [e  e * ft  cos^-e^  + y £ *cos  eA 

- c s it\  9C  hc  e 


E J>  e tA  sen  &D  Mp 

~ CoS  ® S(>»  £' 


l4  f 

£ cos  6 ^ +•  E t> 


.‘4 


CoS  6^  £ CoS 


drt 


Iftt 


+ L 


- 2(Wa  + E £*)/*  + 2si*(Bfi-9^)  B E * 

+ J2  E*  sih  (e„  - <9C)  e"tA  A7C  C 

+ *IV  + E E*  cos  (ep  - 0A)J  e~L*  MpJ> 

+ z[  E VY  Cos  0 + x Sc'h^)  - WAZ  j 


(7) 


Equation  (7)  is  comprised  of  a term,  linear  in  $ with  a known  coef- 
ficient, (a  term  which  would  be  zero  under  either  of  the  two  interference 
hypotheses  mentioned  previously}  two  terms  linear  in  Mc  and  Nhp  with 
coefficients  depending  on  the  unknown  9^  , and  a term  in  which  the 
only  unknown  is  &A  . 


}Z(> 


~ niiMfc 


Similarly, 


rs/ X 

4 t- 

^ bK  ^b^ 


= 2Lc°&*  i-  X cos  i9a  + £ p e~uAsih  6 pc°s  && 


^ £ ft  £it)  <9fl  cos  i9a  +■  £ cos  0g  cos  % e 
+ 2 [ sth * 9 R B - i Sin  8 R - E t>  z'*'* cos  8 -p,  sC*  8 ^ hp 


£ ft  CoS  0A  Sch  (3t>  + C stn  <9j  Sch  -J 


^ 'fcv 


h L— — 

H; 


2 B 4 Ssoh  (0*  - 06)  £ * + -Zcos  £<9B  - /0c)e“‘*A1c  c 

+ 2 (6>»  - &»)  £ + 2 (*  cos  6- S si*  &&) 


Equation  (8)  is  linear  in  £>  , and  A?  , but  with  a coefficient 

which  is  zero  under  the  first  interference  hypothesis.  [Eqn.  (8)  is 
irrelevent  under  the  other  hypothesis].  The  equation  is  linear  in 
and  M p with  coefficients  depending  on  the  unknown  Sg  , and  has  one 
term  in  which  the  only  unknown  ia  9 g . 

Proceeding  to  find  the  derivatives  of  NZ  with  respect  to  9^ 

and  <9B  , 


Mi 

»»« 


aw* 

H 


Similarly, 

bNx 

ie6 


2 £ £ ^/1  /}  * Strt  cos  0fl 


+ 2 He 

[x  F*fi*  t £ fl  B*cosSa  <■  £;>  e-‘4su,es«J)/»Vj'«sl| 

+ 2 Re; 

[ C c.os  £c  e~uA  E*  h ¥ j 

CoS  & /) 

- 2 E 

E * h f)*  s in  8fl  ca  $ 

~ 2 Re 

[ Y E*  /)*  - £ A B * si*  <9fi 

+■  E J>  e-cA  cos  A* E*^Sun <% 

+ 2 Re  ; 

\c  si*  Bc  Mc  e~U  £*/)*/ 

Sch 

2 Re  \ 

A B *E  J cos  (6r  - £„) 

-f  Z Re^fl*  E*  [X  + E J>  e ^ sth  Bp  +■  c e Cos  Bc  ^cj]CoS^ 

+ 2£e.  [fl*  E*  [-Y  ~ E Jp  e cosBphp  + c e^^sio  <9  A?c ]Js^fy 

(9) 


- 2 E>  B*  so*  &B  c°s  &B 

- 2I?€[XB^  + £/9  sih9ft3*  + £ p'e-1-*  s^BpMp  B* 

+ C cos  @c  e V)c^>  * j sui  £?s 
+ 2 3 B*  sth  cos  $g 

- j2fc?  £/  B*t  E /i  c0S  £fl  B*  + £J>  Q ~l'A  cos  6p  Mp  B* 

~ C sih^c  C Mc  B *1  cos 


*9j*wi  winwniiimnwji-inJt,.  • 


= " ^ h 8*)  cos  (0ft  - 0B) 

~ 2 fe£t>*[V  - C sit\  9C  z.  + EP  e ^^cos  6^  Mp  jjcos 

2Re[h  [X  + c cos  &c  e C AMc  + £ J>  e si*  9j>  MjJlJ  sc*  &g 

(10) 


One  hypothesis  models  interference  as  Rayleigh  and  Love  waves 
from  a known  epicenter.  Under  this  hypothesis,  - <9g  = Bx  , and 
Eqs.  (9)  and  (10)  become  irrelevant. 

Alternatively,  if  interference  is  Rayleigh  wave  motion  from  an 
unknown  direction,  Lq.  (10)  becomes  irrelevant,  and  Eq.  (9)  can  be 
simplified  by  using 


B 


* 


o 


(ii) 
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3.  0 Application  to  Rayleigh/ Love  Wave  Interference 

Assuming  that  the  interference  consists  of  a combination  of 
Rayleigh  and  Love  waves  from  a known  direction  B - Q ^ - 6&, 
the  equations  derived  above  may  be  applied  to  find  a minimum- noise 
solution  for  the  signal  magnitude,  and  also  for  the  interfering  wave 
amplitudes.  Equations  (7)  and  (8)  may  be  set  equal  to  zero,  and  by 
introducing  the  subscript  ft,  , for  frequency,  Eq.  (6)  may  be  summed 
over  frequency  and  set  equal  to  zero.  For  simplicity,  «•}  is  symbolic 
of  various  diffe rent  functions  which  may  be  defined  by  reference  to 
Eqs.  (6),  (7),  and  (8),  and  K is  symbolic  of  a subset  of  the  known 
variables.  The  modified  Eqs.  (6),  (7),  and  (8)  are 

[£s„(k)]m  +■  +SjK)r>l  0 (12) 

ft,  h 

St/k^n  + =°  (13) 

Bn  + 5h(K)M  + Sh0<)  = O (14) 

Solving  Eq.  (13)  for  , and  Eq.  (14)  for  Bn  , and  sub- 

stituting the  complex  conjugates  of  these  results  into  Eq.  (12)  gives 

£ZSh(k)  ]m  + 2 Re  [sJk)  f Sh^)w]  <15) 

n.  ^ 

which  is  the  same  as 

[>oSh{«)]^  f ZS„(k)  =0  (16) 

h.  ^ 
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This  symbolic  equation  has  a very  simple,  closed-form  solution 
for  the  minimum-noise  value  of  M . 

Having  solved  for  M , it  is  possible  to  evaluate  further 
minimumization  of  noise  which  might  result  from  small  changes  in 
the  "known"  variables  6 , X , S , and  the  ratio  Results 

corresponding  to  alternate  values  of  the  "known"  variables  may  be 
substituted  into  Eq.  (2)  via  Eq.  (3a- c),  and  those  values  producing 
minimum  T,n%  may  be  taken  as  revised  estimates  of  the  "known 


variables. 


4.  0 Application  to  Rayleigh  Coda  Interference 

Following  a similar  procedure,  the  results  can  be  applied  to 
Rayleigh  coda  interference.  In  this  situation,  the  interference  is 
assumed  to  be  Rayleigh  waves  from  unknown,  frequency-dependent 
directions.  In  this  formulation,  2 may  be  set  equal  to  zero,  and 
Eq.  (6)  may  be  summed  over  frequency.  Eqs.  (6),  (7),  and  (9) 
become,  symbolically, 

[ZSJkAh  i -2Ze{Sn(K)+  5Jk)cos8^  +*„(/<)*'» 

+ Sh  (K^m)  cos&fi^  t tSnCK>i^)‘0 

He  [ft  * [ Sh  ( K}  m)  cos  9 ^ t Sa(k}  M)  sJh  - ° 

While  the  number  of  unknowns  and  the  number  of  equations 
represented  by  Eqs.  (17),  (18),  and  (19)  are  equal,  a closed-form 
solution  is  difficult  to  obtain.  Instead,  using  an  assumed  value  of 
M in  Eq.  (18),  h ^ may  be  solved  for  in  terms  of  6^  ^ and  sub- 
stituted into  Eq.  (19),  which  may  then  be  solved  for  6f\^  • After 
substituting  0^  ^ back  into  Eq.  (18)  to  find  ft  n , Z N ^ can 
be  evaluated  in  Eq.  (2).  This  process  may  be  repeated  with  differ- 
ent values  for  M until  the  minimum  value  of  /Hi  ^ ^ is  found. 

tv 

Hence,  the  minimum-noise  value  of  M can  be  found  by  a 
simple,  single- variable  search. 


I 


A 


APPENDIX  A: 


Carve  Fitting  of  the  VLP  System  Response 


A.  1 Sumina  ry 


The  published  data  on  the  Konigsberg  (KON)  very  long  period 

15 

system  response  was  examined  and  curve  fitted  by  a cascade  of 
three  identically  tuned  circuits.  The  amplitude  vs  period  fit  is  based 
on 


ft  = 


3u 

/a* 


in  which 

"T  is  the  period 

Vo  is  the  resonant  period 

B = 2 ir>  , k>  is  the  damping  constant, 
o co 

As  shown  in  Figures  21,  22,  and  23,  the  resulting  fit  is  quite 
good  in  spite  of  the  fact  that  a complicated  instrument  is  so  simply 
represented. 
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Fig.  21:  Triple  Tuned  Model, 

Curve  Fitting  of  System  Response 
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t-r, 


A.  2 Method  of  Curve  Fitting 


The  response  of  the  seismometer  and  galvanometer  alone  is 
given  by 


Rfa)  = 


- Jh  ** 


(/>*  + 2 k(  LOj  A ■♦-  Wj2  X + 2 t*z  X)  + (J,*) 


in  which  ( A>,  , u>{  ) and  ( , Ui^)  are  the  damping  constants  and 

resonant  frequencies  of  the  seismometers  and  galvanometers.  The 
amplitude  vs  period  response  is 

£2  s K 2 T* 

L(i  - r%*y + * 4 ^ r%/] 

The  parameters  of  the  vertical  seismometers  are 

/t{  & 30 

'i, x * 100 

hi  i & h>2  1 

This  response  is  shown  in  Figure  24. 

Comparison  of  the  seismometer-galvanometer  response  with 
the  overall  response  shows  that  the  degree  of  extra  filtering  is 
considerable.  While  this  is  not  unexpected,  it  does  indicate  that 
one  might  as  well  fit  the  overall  response  as  try  to  imitate  the 
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H3*t? 


multi-stage  filtering  process  with  a mathematical  expansion.  For  this 
reason,  a direct  synthesis  of  the  overall  response  was  made. 

Examination  of  three  axes  shows  that  the  amplitude  curves 
are  characterized  by: 

1.  Geometry  symmetry  in  period, 

2.  Skirt  selectivities  of  6 to  9 dli  per  octave. 

The  simplest  network  (fewest  parameters)  capable  of  fitting  this 
type  of  bandpass  is  three  cascaded  identical  single-tuned  circuits. 

A series,  single-tuned  circuit  has  the  impedance  function 


Z (a>)  = R + aL  4- 


and  an  admittance  function 

'i(xi)  = 


x>c 


A2LC  4 4 I 


^/A 


4 


With 


wo 


o = 

the  form  becomes 


'/JTc 

ft/L 


A C 

the  resonant  frequency 
defining  the  damping  factor 


Y ( />) 


jb2  4 2K0u0  x>  + uj0 


A more  symmetrical  form  is 
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w-v.,  ...  - ..  ■ . 


.* 


V Y >o  ) oL 


(%o  + “%)  + 2*-t 


which  translates  into  period 

. 21 r 


A 


w. 


i T 

27 r 


y(j«)  <* 


-2A.o  7-  j(%0  - “%) 


i 

***  + i (*-/r  - r/e.) 

The  corresponding  amplitude  function,  » for  a single-tuned 

circuit  is 


fl.fr)  = 


n'A 


/X*  7-  (Vr  - r/t,  y] 


For  three  cascaded  single-tuned  circuits,  the  result  is 

fiC?)  ; 


1>/  + lxir  - %JT% 


For  convenience  let 


£ = 2X 


Then, 


fl(T)  ‘ [b*4-  fVr  - Vt.TiVx 
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is  the  amplitude  response  of  three  cascaded  single-tuned  circuits, 
norniEli^ccl  to  ct n of  unity  at  the  peak  ZQ  --  T. 

The  procedure  for  curve  fitting  is: 


Amplitude  vs  Period  is  read  off  the  published  data. 

Geometric  mean  period  is  calculated  for  each 
amplitude  level. 

The  average  geometric  mean  period  is  selected 
as  'tg  . Some  improvement  can  be  made  in  the 
fit,  if  desired,  by  ignoring  the  near  peak  values 
which  are  not  accurate. 

With  known,  a point  estimate  of  S can  be  obtained 
from  each  pair  of  amplitude  and  period  values 


These  are  then  averaged  to  give  the  estimate  for  3 
Some  impro\  ement  is  obtainable  by  ignoring  those 
values  corresponding  to  points  near  the  peak. 

If  desired  a minimum  mean  squared  error  fit  could  be  made 

between  the  data  and  the  response  of  the  three  single-tuned  circuits. 

This  is  unnecessary,  however,  to  prove  that  the  data  can  be  so  fitted. 
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APPENDIX  B:  Leakage  Effects  in  Evaluating  the  Spectrum 

and  Analytic  Signal  of  VLP  Seismic  Events 


B . 1 Summary 


Tlae  evaluation  of  spectrum  and/or  the  analytic  function  of  the 
seismic  signal  is  essential  to  the  operation  of  mixed  event  seismic 
processors  such  as  the  Ellipticity  Filter.  The  difficulties  are 
accentuated  for  VLP  signals  because  the  negative  frequency  region 
can  easily  be  smeared  into  the  positive  frequency  region  to  give 
spurious  results. 

While  the  use  of  a uniform  fader  (window)  is  almost  essential 
in  the  evaluation  of  the  analytic  function,  the  leakage  effects  are 
dramatic  unless  very  long  apertures  are  used.  On  the  other  hand, 
when  only  spectrum  is  desired,  the  cosine  fader  can  be  used  and 
this  weighting  function  permits  small  apertures  to  be  used. 
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B.  2 Di 


scussion 


One  of  the  main  difficulties  to  avoid  in  processing  the  very 
long  period  seismic  signals  is  leakage  which  arises  from  the  trun- 
cation effects  of  the  Dll  process.  This  leakage  acts  not  onl^  to 
smooth  the  spectrum,  as  is  usual  for  all  DFT  processes,  but  be- 
cause the  signal  content  is  of  such  low  frequency,  it  also  smears 
the  negative  frequency  content  into  the  positive  frequency  region. 

The  effect  of  negative  frequency  smearing,  when  a uniform 
fader  is  used,  can  be  determined  from  Figures  25  and  26.  For 
the  system  to  have  a dynamic  range  of  40  dB,  a single  DFT  evalu- 
ation requires  an  fT  product  of  16;  while  a double  DFT  ratio  evaluation 
(such  as  elliptic ity)  requires  an  fT  product  of  32.  The  corresponding 
apertures  for  a frequence  component  whose  period  is  as  large  as  60 
seconds  are  960  and  1820  seconds.  These  apertures  are  quite  long 
and  the  signal  multipath  cannot  be  expected  to  be  stable  in  azimuth 
or  in  content  over  this  duration. 

The  Hanning  or  Cosine  Fader,  on  the  other  hand,  as  can  be 
seen  from  Figures  25  and  26,  require  an  fT  product  of  1.  5 and 
1.97  respectively.  The  corresponding  apertures  arc  then  only  100 
and  120  seconds  respectively  which  is  quite  reasonable. 

The  first  conclusion,  therefore,  is  that  the  Hanning  or  Cosine 
Fader  should  be  used  for  spectrum  evaluations  and  thus  provide  the 
inputs  for  the  frequency  domain  realizations  of  the  Ellipticity  Filter. 


1^3 


fT  Product 


Note,  however,  that  the  key  advantage  of  the  uniform  fader 
is  that  the  direct  and  inverse  operations  in  sequence  yield  the  original 
input.  This  is  not  true  for  the  Cosine  Fader.  Thus,  if  one  is 
instrumentizing  time  domain  realizations  of  the  Elliptic ity  Filter 
(using  the  DFT  and  the  anlaytic  function)^  the  Uniform  Fader  is 
preferred.  The  sequence  of  operations  is  then: 

1.  Uniform  Fader  to  obtain  spectrum. 

Apertures  of  1000  seconds  are  required. 

2.  Suppress  the  negative  frequency  output  and  take 
the  inverse  DFT  (frequency  filter  here  if  desired.) 

3.  This  gives  the  analytic  function  as  an  input  to  the 
time  domain  realization  of  the  Ellipticity  Filler. 

4.  If  desired,  the  frequency  domain  realization  can  now 
be  obtained  by  using  reasonably  short  time  apertures 
(200  seconds)  and  using  the  Cosine  Fader  to  reduce 
smearing  between  the  mixed  events.  This  sequence, 
for  the  frequency  domain  realization,  is  cleaner 
(since  there  are  no  negative  frequencies  at  this 
point)  than  using  the  Cosine  Fader  directly,  but  it 

is  probably  not  necessary. 


The  frequency  domain  realizations  can  also  be  done  at  this 
point  using  a short  time  aperture  uniform  fader.  There  will  be 
smearing  between  the  positive  frequency  components  of  the  mixed 
event,  but  the  negative  frequency  smearing  will  be  negligible. 
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B.  3 Selection  of  the  Fader 


The  large  fT  product  required  by  the  uniform  aperture  leads  one 

to  examine  the  benefits  obtainable  by  using  a different  fader  or  aperture 

lb 

weighting  function.  Several  alternatives  exist  (Hannon,  p.273-280) 

The  normalized  windows  are  given  in  Table  V. 

Of  the  faders  given,  the  Tukey-Hanning  has  the  lowest  sidc- 

. 13 

lobes.  This  Hanning  fa.der  (Blackman  and  Tukey,  p.  15  and  p.  98) 
has  the  unnormalized  response 


jut) 


Co;, 


2rrt  \ 
T ) 


K(f ) 


so k rr-FT 
tt-Ft 


j_  sin  7i '(f  - ‘/r)r 

z r[F-  '/t)t 


I T'^rt~r)r 

z r(f^)r 


sin  irfT 

t r f r 


. t rfr  - r 


Sun,  JJ  F~f 

+ 

nfr  + TT  . 


sin,  t Ft  7 r^_ 

t rfr  tr  * - ( rfi 


Thus,  for  large  P 


TT 


(fT)3 


which  is  substantially  better  than  the  uniform  fader.  For  this  reason, 
as  well  as  the  low  sidelobes  present,  the  Hanning  fader  will  be  used. 
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B.  4 


Performance  oX  the  Uniform  Fader 


For  a uniformly  weighted  time  aperatu.re  (fader),  the  leakage 
can  be  evaluated  as  follows.  For  a sinusoid 


XiO  = Cos(oi0t  + 0 ) 


+ + 0"^u;ot  + 


the  analytic  function  is 

x ft)  = e 


+ 


The  effect  of  the  uniform  time  aperature  on  the  analytic  function 


is  to  give  the  spectral  estimate 

x - ±f  e 

+ rk 


j(wot  + -jut 


dt 


e]e' 

~rU 


■ *T 


ctt 


i^o  — u)  yM  - j - lo)  tx 


- e 


r 


i 6° 


jC  U30  - Ul) 

7T  - f)r 


7T  (£  - f) 


-TV-T”1  r'-T.’TT-rc- 


frequency  of  interest  to  be  IIz.  Then, 


X Qu>') 

~ 1 'f-fc 


+ p i Gt  sck,  ZrrfpT 


unfar 


The  effect  is  zero  when  •/JjT  is  integer.  This  is  as  expected. 


For  -fyT  not  integer,  the  maximum  effect  is 


n 


21  ( j 


z-n 


and  the  normalized  result  is  a complex  amplitude  factor  of 


/ ± 


2trf^r 


The  effect  of  this  amplitude  factor  is  phase  angle  dependent.  The 
maximum  amplitude  effect  occurs  for  0^  = 0°, 


Maximum  Amplitude  Error 


-.(it  -1 ) 

\ 2rr&T  ) 


The  maximum  phase  error  occur 


:s  for  & t V; 


% . The  result  ii 


Maximum  Phase  Error 


4-  j l.  -I  I 

- tan  


Two  effects  occur: 


!•  The  signal  itself  is  distorted  in  amplitude  and  phase. 

For  small  errors,  the  maximum  amplitude  distortion  in  dB  is 

Maximum  Amplitude  Error  (dB)  - •>»  , n ( , t-  L ^ 

zwf,r) 

20  lcS‘°  ej^0  ± J~e^) 

^ <oa>o  B 

Zvf,T 

The  maximum  phase  error  in  degrees  is 


Maximum  Phase  Error 


zw-f^r 


radians 


(znY  ffi  T 


degrees 


degrees 


Comparison  of  these  two  results  show  that  a ten  degret 
phase  error  criteria  is  equivalent  to  a 1.5  dB  amplitude  error. 
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2.  The  dynamic  range  of  the  system  becomes  limited  because 
the  larger  signals  obscure  the  smaller  ones.  The  limit 
on  the  dynamic  range  (dB)  is  about 


ZO  locj  -fyT) 


These  results  are  shown  in  Figures  25  and  26. 


J 

1 


5 


I 1 


B.  5 Performanc  e of  the  Manning  or  Cosine  Fader 


For  a sinusoidal  input 


CoS  C 


the  Hanning  Fader  smears  the  negative  frequencies  onto  the  positive 
frequencies  so  that  the  result  is 


X - 


situ  TTffQ-f)  T 


2*  -f-  A 


i-  (fe-f) r 


+ e 


-j  0 iitu  t r(f0tf)r 

> A ° 


it  + nr  i-fr-hn^r* 


Then,  the  amplitude  distortion  factor  at  f - f is 


+ e 


~i20°  sin.  2it-F0T 


Zn-f0  r l-(2f0T): 


Then, 


r 1 

Maximum  Amplitude  Error  = / +•  / 


while 


Maximum  Phase  Error  = i tatu 


Zirf.rO  - 4-f*  T*) 
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As  for  any  fader  there  arc  two  ef  eels,  distinction.  <.  f the  signal 
itself  and  limitations  on  the  dynamic  ran^e. 

1.  Distortion,  of  the  signal  is  evaluated,  as  before,  for  small 
errors. 

Maximum  Amplitude  Factor  5%  2o  /on  ( I ±.  — — ^ 

• **  o7T  ( J 


in  dQ 


2o  loiu>  g 

8 it  c^r)5 

0-SV6 

(f0f)3 


Maximum  Phase  Factor 


9i r(L'T)2, 


radians 


^(f0ry 

7.  /f 

C SoT)f 


degrees 


degrees 


Comparison  of  these  results  shove  that  a ten  degree  phase  error 
is  equivalent  to  about  a 0.5  dB  amplitude  er  >r. 

2.  Dynamic  range  limitations  can  be  obtained  by  noting  that  the 
limit  on  the  dynamic  range  in  dB  is 


Zo  \o^  STT  C-F0T? 
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which  because  of  tha  cubed 
that  for  a uniform  fader. 

These  results  arc  also  shown  in 
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